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Abstract 


The  phenomenon  of  high  speed  impact  is  of  great  interest  to  the  Air  Force  Office 
of  Scientific  Research  and  the  Air  Force  Research  Laboratory’s  Holloman  High  Speed 
Test  Track.  Rocket  sled  tests  at  the  facility  frequently  are  limited  to  velocities  lower  than 
actually  attainable  due  to  damage  to  the  rail  in  the  form  of  gouges.  Direct  observation  of 
this  gouging  phenomenon  is  not  currently  possible.  This  leaves  computational  modeling 
as  the  only  means  to  study  the  phenomenon.  A  computer  model  has  previously  been  used 
to  model  the  development  of  gouging  at  the  Holloman  High  Speed  Test  Track.  However, 
this  model  has  not  been  experimentally  verifiable  due  to  its  complexity. 

This  research  is  primarily  concerned  with  comparing  experiment  and  analysis  of  a 
simplified  gouging  model.  This  simplified  gouging  experiment  utilized  a  30  mm  powder 
gun  to  shoot  cylindrical  projectiles  at  a  target  at  oblique  angles.  Computer  simulations  of 
the  event  overestimated  penetration  depths  by  13  to  29  percent,  which  is  well  within 
acceptable  limits. 

Using  dimensional  analysis,  the  simplified  gouging  model  was  scaled  up  to  an 
equivalent  sled  system  model.  While  this  equivalent  system  does  not  actually  exist,  it 
does  give  reasonable  estimates  for  similar  sized  systems. 
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VALIDATION  OF  A  SCALED  PLANE  STRAIN 
HYPERVELOCITY  GOUGING  MODEL 


I.  Introduction 


Background 

The  United  States  Air  Force  (USAF)  owns  and  operates  a  high  speed  test  track 
facility  at  Holloman  Air  Force  Base  in  New  Mexico  where  they  can  test  various  hardware 
systems  in  simulated  free  flight  conditions.  The  system  consists  of  a  rocket  powered  sled 
that  glides  on  railroad  track  like  rails  as  shown  in  Figure  1.  The  high  velocity  track  has  a 
narrow  gauge  and  is  just  shy  of  ten  miles  long.  The  Holloman  High  Speed  Test  Track 
(HHSTT)  is  operated  by  the  846th  Test  Squadron.  In  2003,  they  set  a  world  land  speed 
record  of  2884.9  m/s.  They  would  like  to  increase  the  speed  to  above  Mach  10  (-3000 
m/s).  In  the  past,  the  HHSTT  has  experienced  serious  problems  with  rail  damage  from  a 
test  run.  The  846th  Test  Squadron  and  the  Air  Force  Institute  of  Technology  (AFIT)  are 
investigating  methods  of  understanding  this  phenomenon  known  as  gouging  in  order  to 
mitigate  the  damage. 
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Figure  1  -  HHSTT  Rocket  Sled 


Gouging  Phenomenon 

The  rocket  sled  is  attached  to  the  rails  via  four  slippers.  The  rail/slipper 
configuration  is  shown  in  Figure  2.  The  rail  material  is  1080  steel,  while  the  slipper  is 
VascoMax  300  (a  high  strength  steel). 


SUPPER  GAP 


Figure  2  -  Rail/Slipper  Configuration 
Gouging  at  the  HHSTT  typically  occurs  at  velocities  greater  than  1.5  km/s. 
Hooser  [1]  found  the  vertical  velocity  to  be  on  the  order  of  1-2  m/s  with  a  1.5  km/s 
horizontal  velocity.  This  shows  that  the  slipper  impinges  the  rail  with  an  angle  of 
approximately  0.03°.  The  resulting  damage  can  range  from  minor  gouges  as  shown  in 


2 


Figure  3  to  total  structural  failure.  The  test  sled  slippers  disintegrate  upon  impact  at  the 
end  of  the  track  and  are  not  recoverable  for  analysis. 


Figure  3  -  Rail  Gouge 


Laird  [2;  3]  was  one  of  the  first  to  investigate  and  define  gouging  for  the  HHSTT. 
Gouging  occurs  when  inertial  forces  are  so  great  that  the  materials  exhibit  fluid  like 
behavior.  The  relative  motion  of  the  bodies  deforms  the  material  that  continuously 
impinges  on  each  other.  The  resulting  plastic  flow  takes  the  form  of  material  jets  which 
continually  grow  and  penetrate  further  into  the  two  surfaces,  initiating  a  gouge  as  shown 
illustratively  in  Figure  4. 
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Figure  4  -  Gouge  Illustration  [2] 

Problem  Statement  and  Objectives 

Laird  and  Szmerekovsky  [4;  5]  studied  the  HHSTT  slipper  and  rail  interactions 
analytically  using  CTH,  a  hydrocode  developed  by  Sandia  National  Labs.  Szmerekovsky 
used  the  actual  test  sled  conditions  to  perform  his  analysis.  Because  of  this,  his  model  is 
accepted  as  the  standard  simulation  of  the  gouging  phenomenon.  Rickerd  [6]  developed 
a  simplified  gouging  model  equivalent  to  the  Szmerekovsky  model  and  simulated  the 
gouging  phenomenon  with  an  oblique  ballistic  impact  using  CTH.  Unfortunately,  this 
simplified  model  requires  velocities  currently  beyond  the  capabilities  of  test  facilities 
here  at  AFRL.  This  research  attempts  to  compare  experiment  and  analysis  of  a  gouging 
model  that  is  possible  here  at  AFIT.  This  model  will  then  be  scaled  up  to  an  equivalent 
HHSTT  sled  system. 
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II.  Theory 


Conservation  Equations 

The  conservation  equations  are  equations  that  express  of  the  laws  of  physics. 
These  equations  are  useful  for  problems  involving  time  dependant  high  loading.  In 
particular,  when  velocities  or  deformations  of  the  structure  are  large.  Virtually  all  fields 
of  mechanics  and  dynamics  are  based  on  these  fundamental  principals: 

(i)  Mass  is  neither  created  nor  destroyed:  "conservation  of  mass" 

(ii)  Rate  of  change  of  momentum  =  Net  force:  "conservation  of  momentum" 

(iii)  Energy  is  conserved,  though  it  may  change  form:  "conservation  of  energy" 
Mathematically,  the  conservation  of  mass  is  stated 

|  pdV  =  constant  (1) 

v 

where  p  is  density  and  V  is  the  volume  of  the  material. 

The  conservation  of  momentum  can  be  stated  in  a  number  of  different  ways. 

Here  mass  is  assumed  to  remain  constant.  This  leads  to  force  equals  mass  times 
acceleration,  or 

17  dV  /0v 

F  =  m —  (2) 

dt 

where  F  is  the  force  applied,  m  is  the  mass  acted  upon,  v  is  the  velocity  of  the  mass,  and  t 
is  the  time  over  which  the  event  occurs.  Another  useful  statement  of  the  conservation  of 
momentum  is  the  impulse-momentum  relation,  which  is  obtained  by  multiplying  both 
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sides  of  the  conservation  of  momentum  equation  by  dt  and  integrating  over  a  period  of 
time,  giving 

/  =  |  Fdt  =  ^mdv  =  mvf  -mv0  (3) 

Where  /  is  the  impulse  applied  over  some  period  of  time  by  the  applied  force,  v/and  v0 
are  the  initial  and  final  velocities  of  the  mass,  and  therefore  the  right  hand  side  of  the 
equation  is  the  momentum  change  over  some  time  period. 

Lastly,  the  conservation  of  energy  equation  written  for  a  discrete  set  of  j  masses  is 

Z  [e0  +\pvI]=Yj{ex  +LpVA  +  W0^  (4) 

j  v  z  J  j  V  z  J 

where  E  is  the  internal  energy  source,  the  14  pv2  terms  are  kinetic  energy.  W()^{  is  the 
work  done  on  the  system  from  the  initial  state,  0,  to  the  final  state,  1  [7]. 

Stress  Waves  in  a  Continuum 

When  pressure  or  load  is  applied  to  a  material,  stress  waves  develop  and 
propagate  within  that  material.  In  continuum  mechanics,  this  fact  is  often  ignored  and  the 
loading  is  said  to  be  either  static  or  quasi-static.  This  practice  works  well  in  most 
structural  analysis  problems  where  pressure  or  loads  are  applied  very  slowly,  but  in  high 
velocity  impact  dynamics  problems,  is  a  poor  simplification  of  what  actually  occurs.  As 
the  affected  particles  accelerate,  they  build  up  compressive  stresses  in  the  particles 
further  from  the  application  point.  The  motion  of  these  particles  is  governed  by  the 
conservation  of  momentum.  The  speed  of  propagation  of  the  stress  wave,  however,  will 
be  shown  to  be  a  material  property. 
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In  the  following  subsections,  three  regimes  of  dynamic  impact  will  be  described 
to  understand  the  effects  of  impact  velocities  on  stress  wave  propagation.  In  the  first 
regime,  elastic  waves  are  formed  by  low  velocity  impacts.  In  the  second  regime  higher 
velocities  cause  both  an  elastic  wave  and  a  plastic  wave  to  form.  The  plastic  wave  trails 
the  elastic  wave  and  will  cause  permanent  plastic  deformation  in  portions  of  the  material. 
In  the  third  regime  the  impact  velocity  is  much  greater  than  the  sound  speed  of  the 
material  which  causes  a  narrow  region  of  discontinuity  in  state  properties. 

Elastic  Stress  Waves. 

It  is  possible  to  calculate  the  elastic  stress  wave  in  a  uniaxial  rod  impact 
experiment.  Figure  5  shows  an  initially  stationary  rod  is  impacted  by  a  rigid,  semi- 
infinite  plate  moving  from  left  to  right  with  constant  velocity  v0,  which  is  less  than  the 
material  sound  speed  c.  After  the  impact,  an  elastic  stress  wave  travels  to  the  right  at  the 
material  sound  speed,  no  matter  what  the  impact  velocity  is. 
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Figure  5  -  Rod  Impact  Experiment 


Behind  the  stress  wave,  location  1,  the  material  particles  have  been  accelerated  to  a 
velocity  equal  to  the  impact  velocity  v0 .  The  momentum  of  that  material  is 
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pAcAtv0 


(5) 


Where  p  is  the  density  of  the  material,  cAt  is  the  length  of  the  rod  which  is  moving  with  a 
particle  velocity  of  v0,  and  A  is  the  cross  sectional  area  of  the  rod.  v0At  is  the  actual 
deformation  imposed. 

With  the  momentum  defined,  an  expression  for  the  impulse  acting  on  the  material 
is  needed  to  define  the  elastic  stress  wave.  From  basic  physics,  the  impulse  is  known  to 
be  the  integral  of  the  force  over  a  period  of  time.  In  the  uniaxial  case,  the  force  is  the 
stress  times  the  area  over  which  it  acts,  giving  an  impulse  of 

(jAdt  (6) 

where  a  is  the  compressive  stress  at  any  point  in  the  bar  occurring  due  to  the  passage  of 
the  stress  wave.  The  elastic  compressive  stress  wave  is  obtained  by  applying  the 
conservation  of  momentum  to  Equations  5  and  6,  assuming  an  infinitesimal  time  step, 
and  dividing  both  sides  by  the  area  and  time  step.  The  result  is 

CT  =  pcVQ  (7) 

where  c  is  the  material  sound  speed. 

The  rod  impact  case  allows  the  material  sound  speed  to  be  determined.  From  the 
rod,  an  infinitesimal  element  of  area  A  experiences  the  passage  of  a  disturbance,  as 
shown  in  Figure  6. 
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dx 

f- - * 


- > 

x 

Figure  6  -  Material  Element 


The  left  and  right  edges  of  the  element  are  at  positions  of  x  and  x  +  dx ,  respectively. 
Assuming  positive  tension  in  the  positive  x  direction,  the  stress  on  the  left  and  right  sides 
of  the  element  are 


O ' %  and 


<T 


da 

x 

dx 


(8) 


respectively.  The  conservation  of  momentum  for  the  infinitesimal  element  is  given  by 


A(-  (Tx)dt  +  A 


cr  +- 


dx 


dt  =  (pA)d 


v 


(9) 


where  the  left  hand  side  of  the  equation  is  the  impulse  applied  due  to  the  stress  on  the  left 
and  right  boundaries,  and  the  right  hand  side  of  the  equation  is  the  momentum  imparted 
over  the  differential  time  step.  Dividing  both  sides  by  Adt  gives 


da 


dv 

=  p — . 
dx  dt 


(10) 


Using  one  dimensional  strain  and  velocity, 


s- 


du 

dx 


and 


V  = 


du 
~dt  ’ 


(11) 
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it  can  be  shown  that 


ds  _  dv 
dt  dx ' 


(12) 


If  stress  is  assumed  to  be  only  a  function  of  strain,  then  Equations  10  and  12  can  be 
combined  to  give  the  uniaxial  wave  equation, 


2  d2u  d2u 
c  — -  =  — - 
dx 2  dt 2 


(13) 


where 


do / 

c2(s)  =  ^i 

p 


(14) 


For  elastic  behavior,  the  numerator  is  the  elastic  modulus.  This  gives  the  bulk  sound 
wave  velocity  for  a  uniaxial  stress  elastic  impact  as, 


C  - 


(15) 


where  Em  is  the  elastic  modulus  and  p  is  the  density.  For  uniaxial  strain  (plane  strain),  the 
bulk  sound  wave  velocity  is 


c  = 


V 


) 


/?(  1  +  f)(1-2f) 


(16) 


where  vis  the  Poisson’s  ratio  [7]  . 

Understanding  how  waves  propagate  through  a  medium  and  what  happens  at  the 
boundaries  will  be  important  in  analyzing  stress  waves  in  subsequent  sections.  Again 
assuming  a  uniaxial  homogeneous  metal  rod  impacting  a  semi-infinite  rigid  wall  at  a 
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velocity  of  v0  perpendicular  to  the  wall  as  in  Figure  7.  Because  the  impact  is 
perpendicular  and  uniaxial,  it  is  assumed  that  there  are  no  three-dimensional  effects.  It  is 
also  assumed  that  the  impact  is  completely  elastic  and  that  the  wave  velocity  within  the 
rod  is  constant.  The  rod  is  moving  to  the  left  at  a  uniform  velocity  of  v0  and  has  zero 
stress  before  impact.  At  the  moment  of  impact  t  =  t0. 

Length  =  L 

v0 

i - 

G=0 

t— to 

Figure  7  -  Rod  Impact  Experiment 

After  impact  (t0  <  t  <  L/c),  a  stress  wave  propagates  to  the  right  at  the  material  sound 
speed  as  shown  in  Figure  8.  Particles  to  the  left  of  the  wave  experience  a  constant 
compressive  stress  as  described  by  Equation  7  and  due  to  continuity  at  the  left  boundary 
(fixed),  the  rod  velocity  must  be  zero.  To  the  right  of  the  stress  wave  the  rod  continues  to 
travel  to  the  left  at  the  initial  velocity  va  because  the  right  side  hasn’t  “felt”  the  impact 
yet.  The  material  sound  speed  limits  how  fast  a  disturbance  can  propagate  throughout  the 
material. 


C 


M  cAt  ► 

r 

v=0 

o=pcv, 

o  o 
>  II 

r 

|  to<t<L/c 

Figure  8  -  Rod  Impact  Before  Reflection 
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At  t  =  L/c,  the  stress  wave  has  reached  the  right  end  of  the  rod,  see  Figure  9.  At  this 
point  the  entire  rod  has  zero  velocity  and  is  under  constant  compressive  stress  defined  by 
Equation  7.  Conservation  of  energy  dictates  that  the  kinetic  energy  of  the  rod  before 
impact  be  converted  into  internal  strain  energy. 


Figure  9  -  Rod  Impact  at  Reflection 


From  any  mechanics  of  materials  textbook  it  is  known  that  a  free  surface,  by 
definition,  cannot  support  an  applied  stress.  The  only  mechanism  that  can  exist  is  a 
reflection  of  the  stress  wave,  due  to  the  free  end.  A  stress  wave  is  set  up,  opposite  in  sign 
but  equal  in  magnitude,  reducing  the  compressive  stress.  After  reflection,  (L/c  <  t  < 
2L/c),  the  right  side  of  the  stress  wave  must  have  zero  stress,  see  Figure  10. 

Conservation  of  energy  dictates  that  the  internal  strain  energy  be  converted  back  to 
kinetic  energy,  and  therefore,  it’s  initial  velocity  v0,  but  to  the  right.  As  the  stress  wave 
reflects  from  the  free  surface,  it  reflects  with  opposite  sign  (tension).  This  tension  wave 
acts  as  an  unloading  wave,  cancelling  the  effects  of  the  incident  compression  wave. 


C 


n 

v=0 

G= pCV 

— »v0 

<7=0 

L/c<t<2L/c 

Figure  10  -  Rod  Impact  After  Reflection 
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At  t  =  2L/c,  the  stress  wave  is  back  at  the  left  end  of  the  rod  attempting  to  apply  a 
tensile  stress  to  the  rigid  wall.  Since  the  tensile  stress  can’t  be  supported  by  the  interface, 
the  rod  then  separates  from  the  wall.  It  rebounds  off  the  wall,  unstressed,  at  the  initial 
velocity  v0,  Figure  11. 

Vo 

- » 

t>2L/c 

Figure  11  -  Rod  Impact  After  Unloading 

Plastic  Stress  Waves. 

For  most  metals,  the  stress-strain  curve  is  characterized  by  a  linear  elastic  region 
at  low  strains  followed  by  a  region  that  may  or  may  not  be  linear  as  well.  When  a  metal 
is  loaded  beyond  the  linear  portion  it  will  become  permanently  deformed.  This  transition 
point  is  known  as  the  elastic  limit,  dynamic  yield  stress,  or  Hugoniot  elastic  limit  [7].  In 
an  impact  experiment  where  the  velocity  is  great  enough  to  generate  a  stress  higher  than 
this  limit,  the  initial  stress  wave  is  an  elastic  stress  wave  which  is  followed  by  a  plastic 
stress  wave. 

There  are  two  basic  theories  for  describing  plastic  waves,  the  rate-independent 
theory  and  the  rate-dependent  theory.  Both  are  uniaxial  stress  and  as  the  names  imply, 
the  theories  differ  in  their  assumption  of  the  importance  of  strain  rate.  The  rate- 
independent  theory  assumes  that  there  is  a  single  stress-strain  curve  that  describes 
material  behavior  and  that  a  material  had  a  bilinear  stress-strain  curve  which  doesn’t 
depend  on  strain  rate.  The  elastic  stress  wave  would  travel  at  the  elastic  sound  speed 
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given  earlier  as,  c  =  — —  ,  and  have  a  magnitude  equal  to  <y=pcv0.  As  stated  earlier,  for 
V  P 


uniaxial  stress,  the  plastic  wave  would  have  a  slower  velocity  of, 


(17) 


where  Ep  is  the  slope  of  the  stress-strain  curve  in  the  plastic  region  and  cp  is  the  plastic 
wave  velocity.  The  magnitude  of  the  stress  wave  would  then  be 

°p  =  Pcpvo  (18) 

where  ap  is  the  plastic  stress  wave  magnitude.  Figure  12  shows  a  bilinear  stress-strain 
curve  and  the  corresponding  wave  profile. 


Figure  12  -  Bilinear  Stress-Strain  Curve  and  Corresponding  Wave  Profile 
An  alternative  analysis  of  the  rate-independent  theory  assumes  that  the  stress- 
strain  curve  is  concave-up  beyond  the  yield  stress  instead  of  linear,  as  in  Figure  13.  This 
generally  occurs  in  states  of  uniaxial  strain  such  as  plate  impact  experiments.  Uniaxial 
stress  stress-strain  curves  are  normally  concave-down. 
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B 


Stress 


0  Strain 

Figure  13  -  Concave-up  Stress-Strain  Curve 
As  the  strain  is  increased  beyond  the  yield  limit,  the  slope  of  the  stress-strain 
curve,  and  therefore  the  velocity  of  the  stress  wave  increases.  This  means  that  the  higher 
stress  waves  will  eventually  overtake  the  low  stress  waves.  At  which  point  a  plastic 
shock  front  is  formed,  as  shown  in  Figure  14  [8]. 


The  rate-dependent  theory  was  developed  to  account  for  strain  rate  dependence  in 
plastic  flow.  This  theory  more  accurately  approximates  the  plastic  stress-strain  region. 
The  most  popular  form  involves  an  overstress  model  and  is  of  the  form 
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(19) 


<7  =  f{s)+  In  1  +  bsp 

where  f(e)  is  the  stress  from  the  quasi-static  stress-strain  curve,  b  is  a  constant,  and  ep  is 
the  plastic  strain  rate.  [8] 

Shock  Waves. 

If  the  impact  velocity  is  much  greater  than  the  material  sound  speed  for  a  uniaxial 
strain  state,  elastic  waves  will  be  overcome  as  discussed  earlier.  This  leads  to  the 
formation  of  shock  waves,  which  are  narrow  regions  in  which  state  properties  vary 
discontinuously.  This  discontinuity  causes  the  conservation  equations  to  break  down 
when  they  are  in  differential  form. 

The  shock  wave  equations  developed  here  are  for  states  of  uniaxial  strain.  The 
equations  are  found  by  applying  the  conservation  equations  to  the  wave.  Take  a 
shockwave,  traveling  from  left  to  right  into  a  semi-infinite  material  that  is  stationary, 
stress  free,  and  has  no  internal  energy,  see  Figure  15.  The  two  states  identified  by 
subscripts  1  and  0  represent  the  physical  state  ahead  of  and  behind  the  shock  respectively. 
The  shock  velocity  is  U,  while  the  local  particle  velocity  is  u. 


Figure  15  -  Conditions  at  a  Shock  Front  [7] 
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First,  the  conservation  of  mass  across  the  shock  states  that  the  mass  flow  entering 
the  shock  must  equal  the  mass  flow  leaving  the  shock.  The  particle  velocity  in  front  of 
the  shock,  u0,  is  zero,  stationary.  Mathematically,  the  conservation  of  mass  is 

podA(uo  -  U)At  +  pxdA(U  -  ux ) A t  =  0  (20) 

where  the  mass  entering  the  shock  from  the  right  is  p0dA(u0-U)At,  and  from  the  left  is 
pidA(U-ui)At.  Eliminating  the  dA  and  At,  we  obtain, 

PoU  =  Pi(U~u  i).  (21) 

Second,  the  conservation  of  momentum  is  developed.  The  change  in  momentum 
across  the  shock  must  equal  the  impulse  applied,  or 

podAUAtux  =  cjxdAAt .  (22) 

This  can  be  simplified  to  give  the  shock  wave  physics  definition  of  stress, 

o-i  =  poUux .  (23) 

Lastly,  conservation  of  energy  says  that  initial  internal  energy  plus  any  work  done 
on  the  system  is  equal  to  the  final  internal  energy, 

IE0  +  =  IEl  (24) 

where  IEo  is  the  initial  internal  energy,  IEi  is  the  final  internal  energy,  and  W0^x  is  the 

work  done  on  the  system  from  state  zero  to  one.  The  internal  energy  is  the  sum  of  the 
internal  energy  source  per  unit  mass,  E,  plus  the  kinetic  energy  of  the  mass,  KE, 

IE 0  =E()  +  KE 0  and  IEX  =E{+  KEX .  (25) 
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The  internal  energy  source  per  unit  mass  can  be  a  chemical  reaction  that  releases  energy, 
or  a  material  with  some  strain  energy  that  is  stored,  or  other  forms.  Initially  the  kinetic 
energy  is  zero,  because  u0  is  zero.  In  the  final  state  the  kinetic  energy  is, 

KE,=^(PoUAt)(u,f  (26) 

where  p„UAt  is  the  mass  of  material  which  is  moving,  and  »/  is  the  velocity  at  which  it 
is  moving.  After  the  shock  has  passed,  the  internal  energy  source  is 

(a  (U-U,  )M)EX  (27) 

where  the  mass  is  found  using  the  velocity  relative  to  the  shock.  The  combination  of 
Equations  26  and  27  yield  the  internal  energy  after  the  shock  has  passed  as, 

IE\  =  ^  (PoUAt)(lli  f  +  (A  (U  -  ux ) A t)Ex  (28) 

The  internal  energy  source  in  front  of  the  shock  is, 

IE0=(poUAt)Eo  (29) 

where  the  term  in  the  parentheses  is  the  mass  of  material  with  the  internal  energy  E0. 
Because  the  mass  in  front  of  the  shock  was  assumed  stationary,  this  also  happens  to  be 
the  internal  energy  of  the  mass. 

As  for  the  work  done  on  the  system,  work  is  known  to  be  a  force  carried  out  over 
some  distance,  which  is, 

W^x=<Jx(uxAt) .  (30) 

If  Equations  28,  29,  and  30  are  combined,  the  conservation  of  energy  equation  for 
a  moving  shock  wave  is 
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(31) 


(poUM)Eo  +  crxuxAt  -  (px  (U  -  ux  )At)Ex  +  i  (j p0UAt){ux  f . 

Eliminating  At  and  rearranging  so  that  the  internal  energy  terms  are  on  the  same  side, 

(A  (U  -  «,  ))£,  -  poUE0  =  or,!/,  -  i (PoU)(ut  )2 . 


(32) 


The  second  term  on  the  left  hand  side  can  be  simplified  using  conservation  of  mass,  and 
the  second  term  on  the  right  hand  side  can  be  simplified  using  conservation  of 
momentum,  giving 


(px (i U  -  ux))Ex  -px(U-  ux )Eo  =  ctxux  - ^ crxux . 

The  right  side  can  be  simplified  and  dividing  through  by  px(U -ux)  which  gives 


(35) 


E)~E0  = 


1/2  <7XUX 
px(U -ux) ' 


(34) 


Then,  solving  the  conservation  of  momentum  equation  for  ut  and  substituting  into  the 
above  equation  gives 


1 

2 


E(~E0=  r 
Px 


ujpxzEA 

Px 


\ 


u- 


Ujfh-Po) 

Px 


\  • 


(35) 


U  can  be  eliminated,  and  finding  a  common  denominator  on  the  bottom  gives 


Ex~E0 


(P\-Po) 


PlPo 


(36) 


Finally,  simplifying  one  step  further  gives  the  commonly  used  conservation  of  energy 
equation  for  shock  waves, 
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(37) 


EX~E0  = 


n 

kPo 


0 

Ay 


This  is  called  the  Hugoniot  relation.  [9] 

The  conservation  equations  by  themselves  present  an  incomplete  picture  of  shock 
waves.  The  three  equations  (Equations  21,  23,  and  37),  which  are  referred  to  as  the 
Hugoniot  equations,  contain  five  unknowns,  U,  u,  E,  a,  and  p.  One  of  the  two  remaining 
independent  variables  is  found  through  the  use  of  a  Hugoniot  curve,  which  is  the  locus  of 
all  attainable  shock  states  that  are  possible  in  a  material.  The  other  independent  variable 
is  found  from  measurement  or  knowledge  of  one  variable  from  the  equation  of  state, 
discussed  in  a  subsequent  section.  A  Hugoniot  curve  is  analogous  to  a  stress-strain  curve 
in  uniaxial  stress,  but  it  is  not  developed  from  one  experiment  that  follows  the  loading 
path  leading  to  equilibrium  from  a  uniaxial  stress  state.  Instead,  the  curve  is  developed 
using  many  planar  impact  experiments,  such  as  Flyer  Plate  experiments  [10],  to  describe 
the  relationship  between  hydrostatic  pressure  and  specific  volume,  see  Figure  16. 


Figure  16-  Hugoniot  Curve  Showing  Loading  and  Unloading  Paths  [7] 


20 


Each  point  in  a  Hugoniot  curve  represents  a  separate  experiment,  and  hence,  each  point  is 
an  equilibrium  state  for  a  specific  experiment.  [9]  In  this  figure,  the  Hugoniot  curve  is 
marked  as  H.  When  an  impact  occurs  the  loading  path  follows  the  line  from  point  A 
where  the  material  has  zero  pressure  but  high  velocity,  to  point  B  along  what  is  called  the 
Rayleigh  line.  The  loading  does  not  follow  the  Hugoniot,  but  rather  occurs  along  a 
straight  line  connecting  the  initial  state  with  the  peak  pressure  of  impact,  Ph,  which  is  on 
the  Hugoniot  curve.  Unloading  occurs  isentropically  along  the  line  marked  S,  which  is 
not  along  the  Hugoniot  curve.  [9] 

Hugoniot  curves  are  only  valid  in  certain  situations  because  they  are  developed 
under  uniaxial  strain  shock  wave  conditions.  However,  another  use  of  Hugoniot  data  is 
that  it  assists  in  the  development  of  equations  of  state,  which  are  more  general.  [7] 
Computer  codes  that  solve  impact  problems  use  an  equation  of  state  to  relate  internal 
energy,  pressure,  and  volume.  Like  a  Hugoniot,  an  equation  of  state  is  developed  using 
planar  impact  experiments. 

Constitutive  Equation 

A  constitutive  equation  is  the  relationship  between  stress  and  strain  used  by 
continuum  mechanics  codes.  As  stated  before,  the  stress-strain  relation  in  most  metals 
produces  two  distinct  areas  of  response  to  a  loading,  quasi-static  and  dynamic.  In  most 
finite  element  codes,  stress  is  assumed  to  be  quasi-static,  which  means  that  the  loading  is 
applied  so  slowly  that  there  aren’t  any  dynamic  loading  effects.  Hooke’s  Law  is  a 
common  quasi-static  model  used  to  understand  the  stress-strain  curve, 

a  =  f(s,EJ  (38) 


21 


where  a  is  the  stress,  s  is  the  strain,  and  Em  is  the  modulus  of  elasticity.  However, 
Hooke’s  Law  will  only  be  used  in  cases  where  the  stress  is  below  the  yield  stress  of  the 
material.  In  cases  where  the  applied  stress  is  greater  than  the  yield  stress,  it  is  necessary 
to  account  for  dynamic  loading  effects.  The  most  common  way  to  account  for  dynamics 
is  to  include  strain  rate  as  a  variable  in  the  constitutive  equation.  In  general  this  becomes, 


a  =  f(s,s,En 


A 


(39) 


where  s  is  the  strain  rate  applied.  In  some  cases,  constitutive  equations  will  also  be  a 
function  of  internal  energy  and  damage.  [1 1]  In  dynamic  problems  under  large  pressure, 
heating  due  to  plasticity  must  be  accounted  for.  Therefore,  temperature  is  required  in  a 
constitutive  equation. 

Constitutive  equations  exist  for  metals,  ceramics,  concrete,  soil,  and  others.  There 
are  also  numerous  constitutive  equations  with  which  stress-strain  behavior  can  be 
modeled.  Only  one  of  which  will  be  used  in  this  study,  the  Johnson-Cook  equation. 

The  Johnson-Cook  Strength  Model  presents  the  von  Mises  flow  stress  as 


f 


g  =  (A  +  Ben)  1  +  Cln£  n-r*w) 


A 


V 


J 


(40) 


where  a  is  the  von  Mises  flow  stress,  s  is  the  equivalent  plastic  strain,  sP  is  the  plastic 

strain  rate  normalized  by  a  strain  rate  of  1.0s'1,  A,  B,  C,  m,  and  n  are  known  as  the 
Johnson-Cook  coefficients.  The  model  accounts  for  temperature  via  the  homologous 
temperature  T*  as 
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rj-i  * 


T  -  T 


T  -T 

melt  room 


(41) 


where  T  is  the  absolute  temperature,  Tmom  is  the  ambient  temperature,  and  Tme/h  is  the 
melting  temperature  of  the  material.  [7] 

There  is  one  minor  disadvantage  to  the  Johnson-Cook  model.  The  model  presents 
strain  rate  sensitivity  as  being  independent  of  temperature,  which  generally  is  not  the 
case.  The  upside  to  this  is  that  by  keeping  strain,  strain  rate,  and  temperature  uncoupled, 
it  is  relatively  straightforward  to  determine  the  coefficients  from  a  few  experiments  at 
various  temperatures  and  strain  rates.  This  fact  will  be  utilized  in  the  subsequent 
sections. 


Equation  of  State 

The  equation  of  state  (EOS)  describes  the  behavior  of  hydrostatic  components  of 
stress  and  strain.  It  can  be  thought  of  as  a  three-dimensional  constitutive  equation  which 
expresses  the  possible  states  a  material  can  achieve.  Equations  of  state  are  needed  to 
model  how  pressure,  density,  and  energy  relate  when  compressibility  effects  and 
irreversible  processes  such  as  shock  waves  are  included  in  the  problem.  In  fact,  when  the 
pressure  is  very  high,  the  EOS  takes  over  for  the  constitutive  equation.  Energy  needs  to 
be  considered  in  the  analysis  because  shock  wave  conditions  result  from  high  strain  rate 
deformation.  [12] 

Stress  and  strain  can  be  broken  down  into  two  components,  the  hydrostatic  or 
volumetric  stress  or  strain  and  the  deviatoric  stress  or  strain, 

H=h]+h]  (42> 
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where  [a]  is  the  stress  tensor,  [ah]  is  the  hydrostatic  stress  tensor,  and  [cd]  is  the 
deviatoric  stress  tensor.  The  hydrostatic  stress,  or  volumetric  stress,  is  associated  with 
volume  changes,  whereas,  deviatoric  stress  is  associated  with  shape  change.  The  two  are 
handled  by  separate  relationships  in  impact  problems.  The  deviatoric  stress  relationships 
were  discussed  in  the  constitutive  equation  section.  The  hydrostatic  stress  relationship  is 
the  EOS. 

The  equation  of  state  can  be  shown  in  a  general  form  of 

E  -  E(P,V)  (43) 

where  E  is  the  internal  energy,  P  is  the  pressure,  and  V  is  the  specific  volume.  An 
alternative  form,  often  used  in  computer  codes  is, 

P  =  P{p,E).  (44) 

There  are  many  different  forms  of  equations  of  states  that  are  commonly  used  in 
solving  impact  problems,  one  of  which  will  be  discussed  here.  The  Mie-Griineisen  form 
is  simple  and  good  for  modeling  high-pressure  shock  related  events.  [10]  It  is  based 
upon  statistical  mechanics,  using  the  energy  of  individual  atoms  to  arrive  at 
thermodynamic  equations.  As  stated  earlier,  the  Hugoniot  is  one  of  the  pieces  of 
information  used  to  develop  equations  of  state.  Here,  the  Hugoniot  pressure  is  used  as  a 
baseline  in  the  Mie-Griineisen  EOS  and  is  given  by  a  third  order  polynomial, 

PH  =Clp  +  C2p2  +C3p3  (45) 

where  Ph  is  the  Hugoniot  pressure,  the  C,  ’s  are  constants,  and  //  is  defined  as 

=  —  (46) 

Po 
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The  C  parameters  are  only  for  cases  where  density  increases.  If  density  decreases,  C2  and 


where  E  is  the  internal  energy  per  unit  mass,  E0  is  the  internal  energy  per  unit  mass  at 
ambient  conditions,  and  T  is  a  constant  called  the  Griineisen  parameter.  This  parameter 
is  assumed  only  as  a  function  of  specific  volume,  and  is  represented  as 


r 


(48) 


In  this  investigation,  CTH  will  be  used  to  analytically  model  gouging.  CTH 
doesn’t  use  an  equation  for  the  EOS.  Instead,  it  uses  a  table  of  pressure,  energy,  and 
density  at  various  states,  called  the  SESAME  model.  There  are  two  advantages  for  using 
a  table,  first,  there  is  no  need  to  calculate  equation  of  state  variables,  and  second,  the 
exact  physical  state  is  used  as  opposed  to  an  assumed  state,  i.e.  a  quadratic  form  as  in  the 
Mie-Griineisen  equation. 


Buckingham  Pi  Theorem 

Dimensional  analysis  will  allow  an  investigation  of  the  local  Rod  Impact  Model, 
both  experimentally  and  analytically  (CTH),  to  be  applied  to  the  sled-slipper-rail  system 
at  the  HHSTT.  The  Buckingham  Pi  Theorem  is  one  such  tool  for  dimensional  analysis, 
and  will  be  used  to  develop  a  set  of  variants  that  can  scale  a  model  so  that  numerical 
results  and  experimental  results  can  be  applied  and  conclusions  drawn  about  the  actual 
physical  structure. 
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According  to  the  Buckingham  Pi  Theorem,  if  a  physical  law  consists  of  a  number 
of  quantities,  {qi}m,  that  have  dimension  and  are  products  and  powers  of  j  independent 
fundamental  dimensions,  Lj,  then  a  unit  free  fundamental  law  can  be  defined  as 

f(qx,q2,qi,---,qm)  =  Q  (49) 


where  m  is  the  number  of  dimensioned  quantities  to  be  used  in  the  analysis  [13].  A 
fundamental  dimension  is  a  quantity  that  is  used  to  describe  a  dimensioned  quantity. 
There  are  many  different  fundamental  systems  that  can  be  used  such  as  FLT  (Force, 
Length,  Time)  and  MLT  (Mass,  Length,  Time).  Take  pressure  for  example,  in  the  FLT 
system,  pressure  would  be  represented  as  FL'2,  in  the  MLT  system,  pressure  is 
represented  as  ML_1T"2.  It  must  be  ensured  that  the  fundamental  dimensions  alone  can 
describe  all  dimensioned  quantities. 

As  mentioned  above,  it  is  possible  to  represent  any  dimensioned  quantity  as  a 
product  of  fundamental  dimensions  raised  to  some  power 


qt  = 


T  d\  T  d2 

1^1  Lj  2 


Ld”  1 

n  J  i 


(50) 


where  qt  is  a  dimensioned  quantity,  Lj  is  a  fundamental  dimension,  and  <4  is  the  power  the 
fundamental  dimension  is  raised  to.  The  dimensioned  quantities  can  then  be  combined  to 
form  invariant  Pi  quantities 


n  =  faf' far  •••(</„)“' 


(51) 


where  the  a-,  ’s  are  a  to  be  determined  exponent.  It  then  follows  that 


Rearranging  this  equation  so  that  all  of  the  L,  quantities  are  together  leads  to 
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where  the  exponents  /?  can  be  described  as 


A 

du 

d\  2 

••  dXm 

ax 

A 

< 

>  = 

d2l 

d22 

d2m 

< 

a2 

> 

A , 

dnX 

dn2  • 

. .  d 

nm  J 

l  m  J 

(53) 


(54) 


Mathematically,  {a}  must  exist  in  the  null  space  of  the  dimension  matrix,  [D]  for  the 
physical  law  to  be  dimensionally  consistent.  This  requires  that  {(3j  =  {()}.  This 
requirement  forces  the  solution  of  Equation  54  to  give  the  products  of  dimensioned 
quantities  that  must  remain  invariant  between  models.  [13] 

Also  according  to  the  theorem,  if  there  are  m  dimensioned  quantities  and  r 
fundamental  dimensions,  then  there  ar  e  k  =  m-r  independent  dimensionless  quantities. 
In  the  MLT  system  there  will  be  r  =  3  fundamental  dimensions.  [13] 


III.  Methodology 

Cinnamon  [14]  concurrently  developed  a  one-dimensional  approach  for  predicting 
penetration  depth  that  is  used  to  predict  gouging  at  the  HHSTT.  He  verified  his  work 
with  a  series  of  four  oblique  ballistic  impact  experiments.  This  work  will  compare 
computer  simulations  of  these  experiments  to  the  measured  depths.  This  is  done  so  as  to 
yield  multiple  studies  from  the  same  experiment.  Then,  through  the  Buckingham  Pi 
Theorem,  this  model  is  scaled  up  to  an  equivalent  HHSTT  sled  system.  While  this 
equivalent  system  might  not  physically  exist,  it  will  give  ballpark  estimates  of  what 
similar  systems  should  experience  under  the  same  conditions. 
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Cinnamon  Experiments 


The  Cinnamon  experiments  consist  of  a  cylindrical  projectile  with  a 
hemispherical  tip  impacting  a  target  at  an  oblique  angle.  The  projectiles  are  fired  from  a 
30  mm  powder  gun  as  shown  in  Figure  17  with  velocities  on  the  order  of  2.2  km/s. 


Figure  17-30  mm  Powder  Gun 

The  projectiles  are  placed  in  a  plastic  sabot  designed  to  separate  midair,  which 
impact  the  sabot  stripper  plate.  The  projectiles  pass  through  the  hole  in  the  center  of  the 
stripper  plate,  enter  the  target  area,  and  impact  the  target,  see  Figure  18.  A  steel  “pusher 
plate”  is  placed  at  the  base  of  the  projectile  to  prevent  the  projectile  from  lodging  into  the 
sabot  during  launch.  Two  high  speed  cameras  digitally  capture  the  event  at  47,000 
frames  per  second. 
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Figure  18  -  Experimental  Target  Area 

Holloman  AFB  supplied  the  projectiles  and  targets,  which  are  an  18  inch  section 
of  rail  from  the  HHSTT.  Figure  19  shows  a  close-up  of  the  target  rail.  The  expected 
direction  of  travel  and  impact  zone  of  the  projectile  is  also  shown.  The  target  orientation 
was  changed  to  match  the  desired  impact  angle. 


Figure  19  -  Experimental  Target 

It  is  worthy  to  mention  that  the  projectiles  fired  during  the  experiment  did  not 
have  desired  flight  stability  as  in  normal  firearms.  Ordinary  firearms  have  a  spiral  grove 
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machined  down  the  bore  of  the  barrel.  The  projectile,  whose  diameter  matches  the 
diameter  of  the  barrel,  will  spin  axially  as  it  travels  down  and  exits  the  barrel.  This 
spinning  stabilizes  the  projectile  during  flight.  The  gun  barrel  used  in  this  analysis  has  no 
grove  and  is  designed  around  the  sabot.  Multiple  sabot  configurations  allows  multiple 
projectile  sizes  to  use  be  fired  from  the  same  gun.  Hence,  the  sabot  and  projectile  did  not 
spin  and  have  a  natural  instability  during  flight. 

Four  tests  were  performed  by  Cinnamon.  Figure  20  through  Figure  21  show  the 
gouge  for  tests  1  and  2.  Direction  of  travel  is  left  to  right  in  all  photographs. 


Figure  20  -  Experimental  Test  1  Gouge 


Figure  21  -  Experimental  Test  2  Gouge 
Tests  #3  and  #4  were  fired  at  a  different  target  rail  of  the  same  material.  Figure 
22  through  Figure  24  show  before  and  after  photos  of  the  target  rail  for  tests  3  and  4. 
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Figure  22  -Target  Rail  for  Tests  3  and  4,  Before  Test 


Figure  23  -  Experimental  Test  3  Gouge 


Figure  24  -  Experimental  Test  4  Gouge 
As  mentioned  earlier,  it  is  assumed  that  the  projectile  is  not  aerodynamically 
stable,  which  leads  to  tumbling.  Figure  21  shows  indirect  evidence  of  this.  Scuff  marks 
are  seen  on  the  upper  side  of  the  gouge  area,  suggesting  that  the  projectile  impacted  at  an 
angle  relative  to  the  axial  centerline.  However,  the  high  speed  cameras  did  reveal  that  the 
projectile  had  a  flat  orientation  relative  to  the  target. 
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It  is  also  worth  mentioning  here  that  the  two  target  rails  have  different  surface 
preparations.  Target  rail  #1  is  uncoated,  while  target  rail  #2  has  an  epoxy  coating.  The 
effect  of  coatings  on  gouging  is  not  within  the  scope  of  this  work. 

Table  1  lists  the  velocities,  impact  angles,  velocity  components,  and  measured 
depths  of  the  experiments.  The  measurements  were  made  by  taking  a  clay  mold  of  the 
gouges.  It  was  attempted  to  “soft  catch”  the  projectiles  after  impact  for  analysis. 
However,  the  projectiles  disintegrated  upon  impact  and  were  not  recoverable. 


Table  1  -  Cinnamon  Experiment  Results 


Parameter 

Test  1 

Test  2 

Test  3 

Test  4 

v  (m/s) 

2225 

2150 

2147 

2163 

Angle  (deg) 

10 

15 

15 

10 

ux  (m/s) 

2191.20 

2076.74 

2073.84 

2130.14 

uy  (m/s) 

-386.37 

-556.46 

-555.68 

-375.60 

Measured  Depth 

0.5 

1.0 

1.0 

0.5 

(mm,  ±  0.1  mm) 

Model  Development 


In  his  work  for  the  HHSTT,  Szmerekovsky  utilized  the  Buckingham  Pi  theorem 
to  scale  his  model  so  that  he  represented  the  actual  physical  problem,  but  on  a  smaller 
scale.  The  Szmerekovsky  Model  is  a  simplified  plane-strain  model  of  the  HHSTT  sled 
system.  It  consists  of  a  typical  sled  system  and  modeled  as  being  evenly  distributed 
across  four  shoes.  The  flow  chart  showing  how  he  simplified  the  system  is  seen  in  Figure 
25. 
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Figure  25  -  HHSTT  Sled  Simplification 

Here  it  is  seen  that  the  model  assumes  %  the  total  mass  on  one  slipper  and  taking  a  unit 
width  slice  to  end  up  with  the  2D  plane  strain  model. 

Szmerekovsky  used  the  conservation  equations  as  the  physical  laws  to  be  applied 
to  the  Buckingham  Pi  Theorem,  as  will  this  work.  Table  2  lists  the  variables  that 
Szmerekovsky  used. 


The  selection  of  variables  is  one  of  the  most  critical  steps  in  applying  the 
Buckingham  Pi  Theorem.  The  more  variables  included,  extraneous  or  not,  the  higher  the 
chances  are  that  some  of  the  resulting  parameters  will  be  costly  or  nearly  impossible  to 
match  between  models.  Most  engineering  problems  involve  certain  simplifying 
assumptions  that  influence  the  selection  of  variables.  A  suitable  balance  between 
simplicity  and  accuracy  is  the  desired  goal.  How  accurate  the  solution  is  depends  on  the 
objective  of  the  study.  Rickerd  performed  a  Buckingham  Pi  analysis  that  resulted  in  nine 
invariants,  two  of  which  proved  impossible  to  match.  This  analysis  will  use  the 
following  reasoning  to  reduce  the  number  of  input  variables,  in  order  to  further  refine  the 
important  parameters  in  hypervelocity  gouging.  If  a  particular  variable  is  a  function  of 
other  variables  that  are  included  in  the  analysis,  then  it  can  be  removed  with  no  decrease 
in  accuracy.  This  is  because  the  Buckingham  Pi  Theorem  will  ensure  proper  scaling  of 


the  component  variables,  and  hence,  the  removed  variable  is  forced  to  be  scaled 
indirectly.  [13]  Table  3  lists  the  changes  to  the  input  variables  for  the  analysis. 


Table  3  -  Buckingham  Pi  Variable  Selection  Changes 


Change 

Variable 

Reason 

Remove 

w 

Plane  Strain  analysis 

Remove 

s 

S  =  0  for  model 

Remove 

E 

E  =  f(m,  ux,  uy) 

Remove 

c 

c  =  f(Em,  p),  p  =  f(m,  1,  h) 

Remove 

P 

p  =  f(m,  1,  h) 

Add 

t 

Include  time  scale 

The  changes  to  the  Buckingham  Pi  input  variables  results  nine  variables  which 
gives  a  total  of  six  invariants  that  must  be  matched  between  models.  The  invariant 
parameter  fl  is  then 


n  =  (mp  (if  (. if  (u,f  (uyf  (<rj‘  (. Emf  (G0f  (tf.  (55, 


34 


Table  4  lists  the  dimensions  of  the  n  parameter  in  terms  of  the  fundamental  dimensions 


M,  L,  and  T. 


Table  4  -  Buckingham  Pi  Dimensioned  Quantities 


Dimensioned  Quantity 

Symbol 

Fundamental 

Dimensions 

Mass 

m 

M 

Height 

d 

L 

Length 

l 

L 

Horizontal  Velocity 

ux 

LT"1 

Vertical  Velocity 

uy 

LT1 

Compressive  Yield  Strength 

® y,c 

ML_1T'2 

Elastic  Modulus 

E„ 

ml't2 

Shear  Modulus 

G0 

ML_1T"2 

Time 

t 

T 

In  fundamental  dimension  form,  Equation  52  becomes 


n=(Mf  (z)“2  (l)“3  (Lr')“4  ( Lrxy  (mtY2)*6  (ML'r2f  (mtT-2)^  (t)*9 .  (56) 


This  reduces  to 


n  =  (m)a  (Lf2  (rf3 

where, 

A  =a1+a6+a7+as,  (57) 

/?2  =  a2  +  a3  +  a4  +  a5  - a6  - a7  - a8, 

J3 3  =  -«4  -a5  - 2 a6  - 2 a7  - 2as  +  a9. 

Setting  the  values  of  /?  to  zero  and  solving  for  m  =  aj,  l  =  a.2,  and  ux  =  0,4  gives 

ax  =  —a6  -a7  -as, 

a2  =  -a2  +  3 a6  +  3  a7  +  3 as  -  a9 ,  (58) 

a4  =  -a5  - 2 a6  - 2 a7  - 2 as  +  a9. 

Rewriting  these  equations  in  vector  form  gives 
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The  columns  of  Equation  59  represent  a  separate  invariant.  The  invariant  is  found  by 
associating  each  dimensioned  quantity  with  its  corresponding  a  value  and  raising  the 
dimensioned  quantity  to  the  power  seen  in  the  column  vector.  In  this  case,  the  invariants 
are  given  by 


The  separate  invariants  are  found  by  setting  one  C  =  1  and  the  others  to  zero,  which  gives 


7V2 

n, 

71, 

X6 


(61) 
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To  maintain  proper  scaling,  these  six  parameters  must  be  matched  in  both  the  Rod 
Impact  Model  and  the  HHSTT  model.  n\  of  Equation  61  defines  the  geometry  aspect 
ratio,  712  defines  the  impact  angle.  713  through  7x5  of  Equation  61  relate  material  properties, 
length,  mass,  and  horizontal  velocity.  7X6  is  a  time  scale  that  can  be  used  to  compare  the 
two  models. 

Rod  Impact  Model. 

The  Rod  Impact  Model  will  now  be  developed,  which  is  based  on  data  from  the 
Cinnamon  experiments.  The  Szmerekovsky  simplified  HHSTT  sled  model,  CTH 
simulation  model,  and  the  Buckingham  Pi  results  are  based  on  a  plane-strain  analysis. 
Theoretically,  a  plane-strain  rectangular  rod  must  be  used  in  the  experiment  for  the 
comparison  to  be  accurate.  One  key  problem,  however,  is  that  most  high  velocity  impact 
experiments  are  performed  using  cylindrical  projectiles,  which  are  axisymmetric.  As  a 
solution,  experimental  results  are  based  on  cylindrical  projectiles  with  a  hemispherical 
end,  axisymmetric,  while  the  CTH  results  will  be  based  on  a  plane-strain  rectangular  rod 
with  a  cylindrical  end.  The  impacting  end  of  the  rod  will  have  a  curved  tip  because  a 
right  circular  cylinder  would  introduce  too  great  of  a  discontinuity  when  the  90°  comer  of 
the  cylinder  struck  the  impact  plate.  Analytically,  CTH  can  perform  an  oblique  impact 
analysis  using  only  a  plane-strain  rectangular  rod.  Therefore,  the  CTH  Rod  Impact 
Model  will  be  based  on  a  plane-strain  rectangular  rod  impacting  a  target  at  an  oblique 
angle.  Figure  26  shows  the  model.  A  subsequent  section  compares  an  axisymmetric  and 
plane-strain  model  to  provide  a  means  of  comparing  experiment  to  computer  simulation. 
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Figure  26  -  Rod  Impact  Model 

The  material  for  the  Rod  Impact  Model  was  chosen  to  be  the  same  as  the  HHSTT, 
forcing  the  material  properties  to  be  constrained.  The  projectile  dimensions,  mass,  and 
material  properties  are  listed  in  Table  5. 


Table  5  -  Rod  Impact  Projectile  Properties 


m 

4.78  g 

i 

2.5  cm 

d 

0.55  cm 

ay,c 

1.447  GPa 

E0 

184.2  GPa 

G0 

71.8  GPa 

Using  the  velocity  components  of  the  Cinnamon  experiments  in  Table  1,  the  Rod 
Impact  II  invariants  of  Equation  61  are  listed  in  Table  6. 


Table  6  -  Rod  Impact  Model  Parameters 


Parameter 

Test  1 

Test  2 

Test  3 

Test  4 

v  (m/s) 

2225 

2150 

2147 

2163 

Angle  (deg) 

10 

15 

15 

10 

ux  (m/s) 

2191.20 

2076.74 

2073.84 

2130.14 

uy  (m/s) 

-386.37 

-556.46 

-555.68 

-375.60 

Tti 

0.220 

0.220 

0.220 

0.220 

-0.1763 

-0.2679 

-0.2679 

-0.1763 

% 

0.9851 

1.0967 

1.0998 

1.0424 

714 

125.4061 

139.6102 

140.0007 

132.6984 

%5 

48.8825 

54.4192 

54.5714 
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Equivalent  HHSTT  Sled  System. 


The  equivalent  HHSTT  sled  system  is  developed  from  the  invariants  of  Equation 
61.  It  should  be  noted  that  the  invariants  by  themselves  do  not  define  a  specifically  sized 
equivalent  system.  Some  information  about  the  typical  HHSTT  sled  system  is  required  to 
enable  scaling  of  the  remaining  components.  For  instance,  the  only  invariants  that  can 
yield  a  scaled  horizontal  velocity  are  712,  ft 3,  and  715.  With  the  material  properties 

constrained,  713, 714,  and  715  govern  the  theoretical  horizontal  velocity  of  the  sled  system. 
With  the  material  property  removed,  the  invariant  becomes 


which  gives  the  horizontal  velocity  as 


Note  that  lowercase  letters  represent  the  rod  impact  model  properties,  while  uppercase 
letters  represent  the  equivalent  HHSTT  system.  The  result  of  Equation  63,  combined 
with  712  of  Equation  61,  will  give  the  theoretical  velocity  components  of  the  equivalent 


system. 

It  is  desired  to  develop  an  equivalent  HHSTT  sled  system  that  is  comparable  in 
size  to  an  actual  sled  system.  Therefore,  L  and  M  in  Equation  63  will  reflect  typical  sled 
quantities.  It  should  be  noted  that  M  is  the  plane  strain  equivalent  mass  as  discussed 
earlier  in  Figure  25,  so  a  typical  slipper  width  is  needed  in  addition  to  L  and  M.  Table  7 
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lists  typical  sled  and  slipper  properties  that  are  used  in  this  work.  The  slipper  height  is 


found  with  the  given  slipper  length  and  tci  of  Equation  61. 

Table  7  -  Typical  HHSTT  Sled  System  Quantities 


Sled  Mass 

200  kg 

Slipper  Length 

20.32  cm 

Slipper  Width 

10.8  cm 

While  it  is  preferable  to  match  a  typical  slipper  height,  the  height  of  the  equivalent 
system  is  governed  by  the  geometry  of  Rod  Impact  model  and  n\  of  Equation  61. 

Figure  27  shows  a  model  comparison.  Both  are  plane-strain  models.  Chapter  IV 
contains  the  complete  model  along  with  the  simulation. 


Equivalent  HHSTT 


Figure  27  -  Rod  Impact  and  Equivalent  HHSTT  Model  Comparison 
As  for  comparing  the  two  models  in  time,  K6  of  Equation  61  governs  the  time 
scale  between  them.  The  time  scaling  is  governed  by 


frcM 

1 1 ) 

l  L  J 

(64) 


again  with  lowercase  being  the  Rod  Impact  model  and  uppercase  being  the  equivalent 
HHSTT  sled  system.  This  gives  the  equivalent  HHSTT  model  time  as 
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(65) 


T  = 


(  l) 

UJ 

W  J 

t. 


CTH  Model  and  Solution  Method 

Solids  that  undergo  high  strain  rates  and  high  energy  impact  can  flow  plastically 
in  some  areas,  yet  still  maintain  deviatoric  stress  properties.  CTH  is  an  advanced 
hydrocode  that  uses  a  finite  volume  algorithm  to  solve  shock  wave  physics  equations  for 
these  kinds  of  problems.  It  was  developed  by  Sandia  National  Labs  in  the  late  1980’s, 
and  was  primary  designed  to  model  3-D,  multi -material,  large  deformation,  shock  wave 
physics  [15]. 

Generally,  there  are  two  different  ways  to  describe  a  continuum.  One  is  a 
Lagrangian  (material)  description,  in  which  the  material  is  divided  up  into  smaller  pieces 
and  the  conservation  equations  are  solved  by  following  pieces  of  the  material.  The  other 
is  an  Eulerian  (spatial)  description,  which  defines  a  volume  in  space  and  solves  the 
continuity  equations  by  tracking  what  goes  through  it. 

CTH  differs  in  the  solution  schemes  from  other  mechanics  codes.  It  employs  a 
two-step  Eulerian  solution  scheme.  The  first  step  starts  with  an  Eulerian  mesh  which  is 
allowed  to  deform  in  a  Lagrangian  manner.  Conservation  of  volume,  mass,  momentum, 
and  energy  are  conserved  in  this  step.  The  governing  equations  are  replaced  with  explicit 
finite  volume  approximations.  The  second  step  remaps  the  deformed  shape  into  the 
original  Eulerian  mesh.  The  volume  flux  between  cells  is  calculated  first  based  on  the 
cell  face  velocities.  A  High  Resolution  Interface  Tracking  (HRIT)  algorithm  is  then 
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implemented  to  ensure  the  right  materials  are  moved  with  the  volume.  Each  material’s 
corresponding  mass,  momentum,  and  energy  is  then  moved. 

CTH  relies  heavily  on  plots  to  present  post-processing  data.  It  has  two  basic  plot 
types  that  will  be  used  in  this  analysis:  material  and  time-history  plots.  Material  plots  are 
one  or  two  dimensional  contour  or  vector  plots  of  various  properties  versus  position  at  a 
single  timestep,  as  seen  earlier  in  Figure  4.  Time-history  plots  are  plots  of  various 
properties  at  predefined  Lagrangian  tracer  points  versus  time.  These  tracer  point 
locations  need  to  be  defined  in  the  pre-processing  phase.  Therefore,  some  knowledge  of 
the  expected  regions  of  interest,  along  with  trial  and  error  is  necessary  for  placement  of 
these  points  and,  ultimately,  meaningful  results. 

CTH  has  numerous  options  and  parameters  at  the  users’  disposal.  Rickerd  [6] 
investigated  the  proper  options  and  parameters  for  modeling  impacts  related  to  the 
hypervelocity  gouging  problem  and  verified  them  experimentally.  These  options  deal 
with  how  to  handle  yield  strength  in  mixed  cells,  material  volume  fractions  and  pressures 
in  mixed  cells,  void  strength,  and  interface  layers.  Rickerd  showed  that  allowing  multiple 
materials  and  pressures  in  a  cell,  with  pressure  relaxation  best  approximates  the  solution. 

CTH  contains  constitutive  models  that  include  strain  and  strain-rate  effects  that 
are  well  suited  for  the  hypervelocity  gouging  problem.  Standard  options  include  rate- 
dependent  models  for  material  strength  formulations,  such  as  Johnson-Cook.  CTH  also 
contains  two  major  EOS  packages  that  are  used  to  investigate  shock-propagation,  among 
others,  which  occurs  in  hypervelocity  impacts.  These  are  the  Analytic  Equation  of  State 
(ANEOS),  and  SNL-SESAME,  a  tabular  EOS. 
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In  this  work,  materials  are  modeled  using  the  SESAME  EOS  package  in  CTH  for 
projectile  and  target  materials  of  VascoMax  300  and  1080  steel,  respectively.  The 
Johnson-Cook  model  for  the  projectile  is  VascoMax  300.  The  Johnson-Cook  model  for 
1080  steel  doesn’t  exist  in  CTH.  As  a  solution,  the  target  is  modeled  as  iron  with  a 
modified  yield  strength  to  simulate  1080  steel.  The  works  of  Szmerekovsky  and  Rickerd 
have  shown  that  this  solution  is  adequate.  Kennan  [16]  and  Yun  [17]  performed  a 
Johnson-Cook  strength  analysis  of  VascoMax  300  and  1080  steels  by  use  of  a  Split 
Hopkinson  Bar  test.  This  data  was  based  on  strain  rates  on  the  order  of  up  to  10  /s. 
These  Johnson-Cook  coefficients  are  used  in  this  analysis  and  are  listed  in  Table  8. 


Table  8  -  Johnson-Cook  Coefficients 


Coefficient 

1080  Steel 

VascoMax  300 

A 

525  MPa 

2170  MPa 

B 

3.59  MPa 

124  MPa 

c 

0.29 

0.0046 

m 

0.7525 

0.95 

n 

0.6677 

0.3737 

CTH  is  capable  of  parallel  data  processing,  which  allows  simulations  to  run  on 
multiple  processors.  All  simulations  in  this  analysis  used  a  cluster  of  64-bit  Linux 
computers,  using  eight  processors.  On  average,  simulations  took  two  hours  to  complete, 
including  post-processing. 


Model  Boundary  Conditions. 

In  all  simulations,  a  semi-infinite  boundary  condition  was  used  along  the 
boundary  of  the  mesh.  When  a  stress  wave  reaches  a  mesh  boundary,  it  is  imagined  to 
continue  unimpeded  and  not  reflect.  Material  boundaries,  such  as  the  projectile  and 
target  edges,  are  treated  realistically,  and  waves  reflect  as  dictated  by  theory. 
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CTH  has  three  options  for  the  treatment  of  material  interfaces,  and  are 
summarized  here: 

1 .  No  Slide  Line:  materials  are  joined  at  the  contact  surface 

2.  Slide  Line:  a  frictionless  sliding  interface  is  implemented  at  the  contact  surface 

3.  Boundary  Layer:  friction  is  included  and  the  sliding  interface  is  shifted  into  the 
softer  material 

In  a  related  study  to  the  hypervelocity  gouging  phenomena  at  the  HHSTT,  Nguyen  [18; 
19]  investigated  the  proper  option,  and  concluded  that  the  No  Slide  Line  option  should  be 
implemented,  and  is  used  for  all  simulations. 

Axisymmetric  -  Plane  Strain  Comparison. 

As  discussed  earlier,  the  Cinnamon  experiments  used  an  axisymmetric  cylindrical 
rod  with  a  hemispherical  end,  whereas  the  CTH  rod  impact  model  is  a  plane-strain 
rectangular  rod  with  a  semicircular  end.  To  compare  the  two,  it  is  necessary  to 
investigate  how  similar  CTH  simulates  an  axisymmetric  and  plane-strain  impact  of  the 
same  configuration.  CTH  can  only  simulate  both  in  a  perpendicular  impact  as  in  Figure 
28.  Therefore,  an  axisymmetric  and  plane-strain  simulation  of  a  vertical  impact  of  the 
same  length  and  width  as  the  projectiles  will  be  performed  in  Chapter  IV.  The  velocity 
will  match  the  y-component  velocity  of  the  actual  rod  impact  experiments. 
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Figure  28  -  Vertical  Impact  for  Axisymmetric  -  Plane  Strain  Comparison 
The  mesh  size  in  the  simulation  was  0.01  cm  square  in  the  region  of  impact  of  the 
target  and  throughout  the  rod.  In  the  area  away  from  the  impact  of  the  target,  the  mesh 
grows  courser  incrementally.  A  diagram  of  the  mesh  is  shown  in  Figure  29,  which  is  not 
to  scale. 


x=0.275  cm 
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For  the  axisymmetric  simulation,  only  half  of  the  model  is  defined.  A  mirror 
option  is  used  during  plotting  to  view  the  entire  model.  For  the  plane-strain  simulation, 
the  entire  model  is  defined  by  mirroring  the  mesh  about  the  left  vertical  edge.  Appendix 
A  contains  the  CTH  input  deck  for  the  axisymmetric  simulation. 

Rod  Impact  Model. 

For  the  rod  impact  model,  the  mesh  was  0.005  cm  square  in  the  area  of  the 
impact.  Farther  from  the  impact  zone,  the  mesh  grows  coarser  incrementally.  The  mesh 
is  shown  schematically  in  Figure  30,  not  to  scale.  The  rod  is  shaded  dark  gray. 


Figure  30  -  Rod  Impact  Model  Mesh 


Appendix  B  contains  the  CTH  input  deck  for  the  rod  impact  simulation  #1.  All  other  run 
are  the  same,  with  the  exception  of  the  velocity  components. 
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IV.  Analysis  and  Results 


In  order  to  simulate  the  Cinnamon  experiments,  the  axisymmetric  -  plane-strain 
comparison  in  CTH  is  analyzed  first.  The  four  Cinnamon  simulations  follow.  The 
results  are  then  compared  to  the  measured  depths  of  the  Cinnamon  experiments.  Using 
the  Buckingham  Pi  analysis,  Test  #1  simulation  is  then  scaled  up  to  the  equivalent 
HHSTT  sled  system  model,  which  is  simulated  in  CTH  as  well. 


Axisymmetric  -  Plane-Strain  Comparison 

The  goal  of  this  investigation  is  to  determine  a  suitable  time  span  that  a  plane- 
strain  CTH  analysis  can  simulate  an  axisymmetric  oblique  impact.  While  the  governing 
equations  remain  the  same  regardless  of  an  axisymmetric  or  plane-strain  analysis,  the 
forms  of  the  equations  are  different.  On  a  local  spatial  scale,  the  differences  can  be 
dramatic,  especially  towards  the  edge  of  the  projectile.  This  is  due  to  the  fact  that  in  an 
axisymmetric  analysis,  the  stress  must  be  zero  at  the  free  end  and  the  stress  wave  reflects 
with  opposite  sign,  whereas,  in  a  plane-strain  analysis,  the  velocity  is  zero  at  the  fixed 
end  and  the  stress  increases.  However,  of  primary  importance  is  how  similar  the 
simulations  can  predict  penetration  depth  while  maintaining  similar  shape,  which  is  more 
on  the  global  scale.  Also  of  importance  is  how  similar  pressure  profiles  are  between  the 
two.  With  this  said,  plots  of  the  displacement  and  pressure  fields  will  be  compared  side- 
by-side.  The  deciding  criteria  for  determining  any  dissimilarity  is  the  time  at  which  the 
displacement  field  within  the  target  interior  starts  to  differ  between  the  two. 


47 


With  the  Cinnamon  experimental  velocities  as  stated  above,  two  comparisons  are 


investigated  using  vertical  velocities  of  375  and  550  m/s.  First,  the  375  m/s  velocity, 
corresponding  to  Test  #1  and  Test  #4,  is  analyzed.  This  is  followed  by  the  555  m/s 
velocity  analysis,  corresponding  to  Test  #2  and  Test  #3.  Figure  3 1  shows  the  simulation 
at  the  moment  of  impact  (t  =  0  s)  for  both  axisymmetric  and  plane-strain.  It  also  shows 
the  location  of  three  Lagrangian  tracer  points  located  at  (0,  0.01),  (0,  0.275),  and  (0.275, 
0.275)  cm,  respectively. 
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Figure  31  -  Axisymmetric  and  Plane-Strain  Comparison  (t  =  0  s) 

375  m/s  Impact  Velocity. 

For  the  375  m/s  velocity  case,  Figure  32  through  Figure  39  show  a  side-by-side 
comparison  of  the  displacement  field  of  the  axisymmetric  and  plane-strain  simulations  at 
times  of  2.5,  7.5,  10.0,  and  12.0  microseconds.  All  figures  on  the  left  are  axisymmetric, 
and  all  figures  on  the  right  are  plane-strain. 
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Figure  32  -  v  =  375,  axi,  t=2.5  jus 


Figure  33  -  v  =  375,  pi  st,  t=2.5  p is 


Figure  34  -  v  =  375,  axi,  t=7.5  jus 


Figure  35  -  v  =  375,  pi  st,  t=7.5  jus 
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Material  Imp-act 


Material  Impact 


2DC  X  (cm) 

Axisym metric  Impact:  V3G0  cm  103 0,  v^37£  m/s 

Figure  36  -  v  =  375,  axi,  t=10.0  /us 


20K  X  (cm) 

Plane  Strain:  V300  an  1080,  y=375  m/a 
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Figure  37  -v  =  375,  pi  st,  t=10.0  /us 

Material  Impact 


X  (cm) 


Axbymmetrfe  Impact:  VJOO  on  1080,  v=375  m/s 

Figure  38  -  v  =  375,  axi,  t=12.0  ps 


St  (cm) 


Plano  Strain:  V300  an  1080,  v=375  m/s 

Figure  39  -  v  =  375,  pi  st,  t=12.0  ps 


From  the  very  beginning,  it  is  seen  that  the  plane-strain  analysis  predicts  a  larger 
amount  of  target  material  to  be  displaced  out  and  upwards  of  the  impact  zone.  While  this 
is  a  significant  difference,  the  conditions  at  the  surface  of  the  target  away  from  the 
projectile  are  not  important.  What  is  important  is  what  happens  in  the  target  interior. 

As  the  simulation  progresses,  a  key  difference  is  seen  in  the  deformation  of  the 
projectile  tip.  At  10  ps,  the  projectile  tip  has  flattened  out  in  the  axisymmetric 
simulation.  The  tip  of  the  plane-strain  simulation  has  already  flattened  and  started  to 
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form  a  concave  down  shape.  At  12  ps,  the  flattened  tip  of  the  axisymmetric  simulation 
has  increased  its  flatness,  while  the  plane-strain  simulation  progresses  further  into  a 
concave  down  shape.  The  penetration  shape  of  the  plane-strain  simulation  starts  to 
deviate  from  the  axisymmetric  version,  suggesting  that  the  projectile  has  started  to  split  at 
the  tip.  This  is  due  to  the  fact  that  the  axisymmetric  is  3D,  whereas,  the  plane-strain  is 
ID.  Because  of  this,  10  ps  is  chosen  to  be  the  maximum  time  limit.  This  implies  that  the 
plane-strain  simulation  of  Test  #1  and  #4  of  the  Cinnamon  experiments  are  considered 
accurate  up  to  10  ps. 

The  time  history  plots  of  CTH  show  what  happens  on  a  local  spatial  scale 
throughout  the  event.  Figure  40  and  Figure  42  show  y  position  as  a  function  of  time  of 
tracer  points  (TP)  1,  axisymmetric  on  the  left. 


TIMF  Ois) 

Axisymmetric  Impact:  V30U  on  1080,  *=375  m/s 

Figure  40  -TP  1  y  Position ,  axi 


£  «2. 


TWF  (je*s) 

Plain  Strain:  V30D  on  1080,  *=375  m/s 


Figure  41  -  TP  1  y  Position,  pi  st 


Figure  42  and  Figure  44  show  y  position  as  a  function  of  time  of  the  projectile  edge,  TP 
3. 
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Plain  Strain:  VSOO  on  TQBD„  v-375  m/s 

Figure  43  -  TP  3  y  Position,  pi  st 


Figure  42  -  TP  3  y  Position,  axi 
It  is  seen  that  here  that  the  two  simulations  are  similar  up  to  roughly  10  ps  further 
illustrating  that  the  simulations  are  quite  similar  up  to  this  point  in  time. 

With  the  time  span  defined  above,  the  pressure  profiles  are  now  compared. 
Figure  44  through  Figure  49  show  side-by-side  comparison  of  pressure  profiles  at  times 
of  1.5,  5.0,  and  10.0  microseconds.  Again,  axisymmetric  on  the  left. 
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Figure  44  -  Pressure,  axi,  7=7.5  jus 
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Plane  Strain:  V300  on  1080.  v=375  m/s 


Figure  45  -  Pressure,  pi  st,  t=1.5  ps 
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Figure  47  -  Pressure,  pi  st,  t=5.0  jus 


Pressure 
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Asymmetric  Impact:  V300  on  1030,  v=375  m/s 


?0R  X  (cm) 

Plane  Strqfn:  V300  On  108 D,  v=375  m/s 


Figure  48  -  Pressure,  axi,  t=10.0  ps  Figure  49  -  Pressure,  pi  st,  t=10.0  ps 

It  is  seen  that  the  two  simulations  differ  in  the  profile  development,  but  the  relative 
difference  in  magnitude  is  not  too  severe.  To  further  see  the  differences,  Table  9  shows 
the  relative  difference  of  pressure  at  six  points  (in  Eulerian  space).  Two  x  positions  of  0 
and  2.5,  and  three  y  positions  of  -0.2,  -0.4,  and  -0.6  are  chosen  sample  points,  and  shown 
on  the  last  two  figures. 
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Table  9  -  Pressure  Difference,  v  =  375  m/s 


Percent  Difference 

X 

y 

t=1 .5  (J.S 

l-K 

II 

Ol 

■p 

</> 

t=10  jus 

0 

-0.2 

14% 

21% 

16% 

0 

-0.4 

28% 

9% 

14% 

0 

-0.6 

53% 

0% 

16% 

0.25 

-0.2 

52% 

17% 

30% 

0.25 

-0.4 

19% 

47% 

0% 

0.25 

-0.6 

19% 

36% 

12% 

avg 

31% 

21% 

14% 

This  leads  to  the  conclusion  that  pressure  profiles  in  the  Cinnamon  simulations  of  Tests 


#1  and  #4  will  inherently  have  differences  on  average  of  14  to  31  percent. 


555  m/s  Impact  Velocity. 


For  the  555  m/s  velocity,  Figure  50  through  Figure  57  show  a  side-by-side 
comparison  of  axisymmetric  and  plane-strain  simulations  at  times  of  2.0,  6.0,  8.0,  and 
10.0  microseconds.  All  figures  on  the  left  are  axisymmetric,  and  all  figures  on  the  right 


are  plane-strain. 

Material  Impact 
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Axisymmetric  impact:  V3QQ  cm  1030,  V=555  m/s 

Figure  50  -  v  =  555,  axi,  t=2.0  jus 


2  DR  X  (cm) 

Plane  Strain:  V30D  on  10RDT  v=555  m/s 

Figure  51  -v  =  555,  pi  st,  t=2.0  ps 
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Figure  52  -v  =  555,  axi,  t=6.0  jus 


Figure  53  -  v  =  555,  pi  st,  t=6.0  jus 
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Figure  54  -v  =  555,  axi,  t=8.0  jus 


Figure  55  -  v  =  555,  pi  st,  t=8.0  /us 
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Figure  56  -  v  =  555,  axi,  t=10.0  jus  Figure  57  -  v  =  555,  pi  st,  t=10.0  /us 

As  with  the  375  m/s  velocity,  from  the  very  beginning,  it  is  seen  that  the  plane- 
strain  analysis  predicts  a  larger  amount  of  target  material  to  be  displaced  out  and  upwards 
of  the  impact  zone.  Also  as  with  the  375  m/s  velocity,  another  difference  is  the  manner 
in  which  the  top  surface  of  the  target  away  from  the  impact  deforms.  The  ‘humps’  of  the 
plane-strain  simulation  are  more  pronounced  and  spread  out.  While  these  are  significant 
differences,  the  conditions  at  the  surface  of  the  target  away  from  the  projectile  are  not 
important.  What  is  important  is  what  happens  in  the  target  interior. 

As  the  simulations  progress,  a  slight  difference  is  seen  in  the  deformation  of  the 
projectile  along  the  edge,  while  the  tip  remains  circular.  While  the  plane-strain  shape  is 
not  a  curved  as  the  axisymmetric,  it  does  not  worsen  as  time  progresses.  At  8  ps,  the 
projectile  tip  starts  to  flatten  out  in  the  plane-strain  simulation.  At  10  ps,  the 
axisymmetric  simulation  still  maintains  a  flattened  tip,  while  the  plane-strain  simulation 
has  started  to  form  humps  at  the  tip.  The  penetration  shape  of  the  plane-strain  simulation 
starts  to  deviate  from  the  axisymmetric  version.  Because  of  this,  10  ps  is  chosen  to  be 
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the  maximum  time  limit.  This  implies  that  the  plane-strain  simulation  of  Test  #2  and  #3 
of  the  Cinnamon  experiments  are  considered  accurate  up  to  10  ps. 

The  time  history  plots  of  CTH  show  what  happens  on  a  local  spatial  scale 
throughout  the  event.  Figure  58  and  Figure  60  show  the  y  position  of  tracer  point  1, 
axisymmetric  on  the  left. 


TIMF  (fri) 

Ax  Symmetric  Impact:  V3DG  on  1080*  m/s 


TIME 

Plain  Strgln:  V30D  on  1080.  m/s 


Figure  58  -TP  1  y  Position,  axi  Figure  59  -TP  1  y  Position,  pi  st 

It  is  seen  from  these  two  figures  that  the  two  simulations  are  similar  up  to  roughly  12  ps. 
Figure  60  and  Figure  62  show  the  y  position  of  the  projectile  edge,  TP  3. 
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TIME  (/is) 

Am  [symmetric  Impact:  V3D0  on  1080.  v=555  m/s 

Figure  60  -TP  3  y  Position,  axi 


TIMF  (>s) 

Plain  Strafn:  Y300  on  10  BO,  m/s 

Figure  61  -  TP  3  y  Position,  pi  st 


It  is  seen  that  here  that  the  two  simulations  are  similar  up  to  roughly  10  ps  further 


illustrating  that  the  simulations  are  quite  similar  up  to  this  point  in  time. 

With  the  time  span  defined  above,  the  pressure  profiles  are  now  analyzed.  Figure 
62  through  Figure  67  show  side-by-side  comparison  of  pressure  profiles  at  times  of  1.5, 


5.0,  and  10.0  microseconds.  Again,  axisymmetric  on  the  left. 
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Figure  62  -  Pressure,  axi,  t=1.5  ps 


Figure  63  -  Pressure,  pi  st,  t=1.5  ps 
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Figure  64  -  Pressure,  axi,  t=5.0  jus 
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Figure  65  -  Pressure,  pi  st,  t=5. 0  ps 
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Figure  66  -  Pressure,  axi,  t=10.0  ps 


2DR  X  (cm) 

Plane  Strain:  V30D  on  ffiSCT  v=555  rr/s 

Figure  67  -  Pressure,  pi  st,  t=10.0  ps 


It  is  seen  that  the  two  simulations  differ  in  the  profile  development,  but  the 


relative  difference  in  magnitude  is  not  too  severe.  To  further  see  the  differences,  Table 


10  shows  the  relative  difference  of  pressure  at  six  points  (in  Eulerian  space).  Two  x 
positions  of  0  and  2.5,  and  three  y  positions  of  -0.2,  -0.4,  and  -0.6  are  chosen  sample 


points,  and  shown  on  the  last  two  figures. 
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Table  10  -  Pressure  Difference,  v  =  555  m/s 


Percent  Difference  j 

X 

y 

t=5  jus 

0 

-0.2 

74%  1 

0 

-0.4 

0 

-0.6 

10% 

0.25 

-0.2 

16% 

7% 

0% 

0.25 

-0.4 

52% 

42% 

13% 

0.25 

-0.6 

52% 

42% 

0% 

avg 

39% 

29% 

14% 

This  leads  to  the  conclusion  that  pressure  profiles  in  the  Cinnamon  simulations  of  Tests 


#2  and  #3  will  inherently  have  differences  on  average  of  14  to  39  percent. 


Bringing  this  all  together,  Table  1 1  shows  a  time  range  comparison  between  the 


allowed  contact  time  of  the  experiment  and  the  time  limit  imposed  on  the  CTH  analyses 


from  above.  The  contact  time  was  taken  to  be  the  time  to  transverse  one  full  length  of  the 


projectile  traveling  at  the  horizontal  velocity. 


Table  1 1  -  Time  Range  Comparison 


Test 

Velocity 

(m/s) 

Angle 

(deg) 

Experimental 
Contact  Time 

/list,  approx. 

Time  of  CTH 
Simulation  (ps) 

Percent 

Simulated 

1 

2225 

10 

10 

88% 

2150 

15 

10 

83% 

2147 

15 

10 

83% 

M 

2163 

10 

WMBEMM 

10 

85% 

It  is  seen  that  CTH  is  capable  of  capturing  most  of  the  event,  but  not  its  entirety. 


Rod  Impact  Simulation 

With  the  imposed  CTH  time  limit  above,  the  simulations  of  the  Cinnamon 
experiments  are  made.  Figure  68  shows  the  configuration  at  the  moment  of  impact  (t  =  0 
5),  for  all  simulations.  As  mentioned  in  the  discussion  of  the  Cinnamon  experiments,  the 
projectile  is  oriented  level  against  the  target  with  the  velocity  vector  at  the  appropriate 
angle,  as  shown.  The  only  difference  between  simulations  is  the  velocity  vector.  The 
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CTH  simulations  of  Tests  #1  and  #4  are  analyzed  first  since  they  share  the  same  impact 
angle  of  10°.  The  analysis  of  Tests  #2  and  #3,  which  share  the  angle  of  15°,  follows. 
The  simulations  are  stopped  at  10  ps  because  the  plane-strain  simulation  of  the 
axisymmetric  event  is  not  considered  accurate  past  this  point. 


Material  Impact 


0.0  0.5  1.0  1.5  2.0  2.5  5.0 

2DG  X  (cm) 


Figure  68  -  CTH  Rod  Impact  Model,  (t  =  0  s) 

Test  #1  Simulation,  v=2225  m/s  @  10°. 


Figure  69  through  Figure  72  show  the  simulation  at  interesting  times  of  3.0,  4.5, 
6.5,  and  10.0  microseconds. 
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Figure  69  -  Test  #1  Sim,  t=3. 0  ps 


Figure  70-  Test  #1  Sim,  t=4.5  ps 
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Figure  71  -  Test  #1  Sim,  t=6.5  jus  Figure  72  -  Test  #1  Sim,  t 

From  3  to  4.5  pis,  the  projectile  and  target  start  do  deform  and  form  a  gouge.  At  6.5  ps, 


the  projectile  is  on  the  verge  of  forming  a  crack.  At  10  ps,  the  projectile  crack  is 

developing,  and  the  penetration  depth  of  the  target  is  roughly  1  cm.  To  determine  a  better 

estimate,  Figure  73  shows  a  blown  up  image  of  the  simulation  at  10  ps.  Note  that  the  x 

and  y  axis  scaling  are  different. 
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Figure  73  -  Test  #1  Sim,  t=10  /us,  zoom 
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It  is  seen  that  the  simulated  penetration  depth  is  on  the  order  of  0.7  cm.  The  measured 
depth  is  0.5  ±0.1  mm.  The  difference  in  depth  is  roughly  29  percent,  which  is  very  good. 
It  is  also  interesting  to  note  pressure  and  plastic  strain-rate  profiles  in  Figure  74 


and  Figure  75.  Due  to  less  difference  in  pressure  magnitudes  at  10  ps,  this  time  is  chosen 


to  better  represent  the  pressure  distribution. 
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Figure  74  -  Test  #7  Pressure  Profile,  t=10.0  jus  Figure  75  -  Test  #7  Strain  Rate,  t=10.0  jus 


It  is  seen  here  that  pressures  on  the  order  of  3.0  GPa  are  present,  indicating  that  the 


plastic  limit  is  exceeded  and  shock  waves  occur.  It  is  also  seen  that  strain-rates  as  high  as 
104  /s  are  present,  indicating  that  although  the  current  Johnson-Cook  strength  model  is 


fairly  accurate,  it  has  just  exceeded  its  range  of  applicability.  As  noted  earlier,  the 
Johnson-Cook  strength  model  was  based  on  strain  rates  up  to  103  /s. 


Test  #4  Simulation,  v=2163  m/s  @  10°. 

Figure  76  through  Figure  79  show  the  simulation  at  interesting  times  of  3.0,  4.5, 
6.5,  and  10.0  microseconds. 
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Materiel  impact 
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Figure  76  -  Test  #4  Sim,  t=3. 0  /us 
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Figure  77  -  Test  #4  Sim,  t=4.5  jus 
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Figure  78  -  Test  #4  Sim,  t=6.5  jus  Figure  79-  Test  #4  Sim,  t=10.0  /is 

From  3  to  4.5  ps,  the  projectile  and  target  start  do  deform  and  form  a  gouge.  At  6.5  ps, 
the  target  material  has  formed  the  classic  jetting  of  material  more  typical  of  HHSTT 
gouges.  At  10  ps,  the  projectile  is  developing  a  crack  on  the  top  side,  and  the  penetration 
depth  of  the  target  is  roughly  0.1  cm.  To  determine  a  better  estimate,  Figure  80  shows  a 


blown  up  image  of  the  simulation  at  10  ps.  Note  that  the  x  and  y  axis  scaling  are 
different. 
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Figure  80  -  Test  #4  Sim,  t=10.0  jus 

It  is  seen  that  the  simulated  penetration  depth  is  on  the  order  of  0.7  mm.  The  measured 
depth  is  0.5  ±0.1  mm.  The  difference  in  depth  is  roughly  29  percent,  which  again  is 
quite  good.  Pressure  and  plastic  strain-rate  profiles  are  similar  to  Test  #1  with  pressures 
strain  rates  up  to  3  GPa  and  104  /s,  respectively,  and  are  shown  in  Appendix  C. 

Compared  to  Test  #1,  a  couple  noticeable  differences  in  the  simulation  are  the 
target  gouge  formation  and  the  projectile  crack.  The  gouge  formation  in  Test  #4  is  more 
pronounced  than  in  Test  #1,  whereas  the  projectile  crack  is  more  pronounced  in  Test  #4 
than  Test  #1.  This  is  due  to  the  fact  that  the  total  velocity,  and  hence  energy,  is  lower  in 
Test  #4  than  Test  #1. 


Test  #2  Simulation,  v=2150  m/s  @  15°. 

Figure  81  through  Figure  83  show  the  simulation  at  interesting  times  of  2.5,  5.5, 
and  10.0  microseconds. 
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Figure  81  -  Test  #2  Sim,  t=2.5  jus  Figure  82  -  Test  #2  Sim,  t=5.5  jus 
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Figure  83  -  Test  #2  Sim,  t=10.0  jus 

From  2.5  to  5.5  ps,  the  projectile  and  target  start  do  deform  and  form  a  gouge.  This 
gouge  is  more  pronounced  than  Tests  #1  and  #4.  Also  at  5.5  ps,  the  projectile  has  formed 


a  crack  on  the  under  surface  just  ahead  of  the  gouge.  At  10.0  ps,  the  projectile  has 
formed  another  crack,  and  both  almost  propagated  all  the  way  through.  The  penetration 
depth  of  the  target  is  roughly  1  cm.  To  determine  a  better  estimate,  Figure  84  shows  a 
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blown  up  image  of  the  simulation  at  10  ps.  Note  that  the  x  and  y  axis  scaling  are 
different. 
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Figure  84-  Test  #2  Sim,  t=10.0  jus 


It  is  seen  that  the  simulated  penetration  depth  is  on  the  order  of  1 .25  cm.  The  measured 


depth  is  1.0  ±  0.1  mm.  The  difference  in  depth  is  roughly  20  percent,  which  is  excellent. 


It  is  also  interesting  to  note  pressure  and  plastic  strain-rate  profiles  in  Figure  85 


and  Figure  86.  Due  to  less  difference  in  pressure  magnitudes  at  10  ps,  this  time  is  chosen 


to  better  represent  the  pressure  distribution. 
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Figure  85  —  Test  #2  Pressure  Profile,  t=10.0  jus  Figure  86-  Test  #2  Strain  Rate,  t=10.0  /js 
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It  is  seen  here  that  pressures  on  the  order  of  3.0  GPa  and  strain-rates  as  high  as  104  /s  are 
present.  As  with  Test  #1  and  Test#4,  this  indicates  that  the  plastic  limit  is  exceeded  and 
shock  waves  occur  and  that  although  the  current  Johnson-Cook  strength  model  is  fairly 
accurate,  it  has  just  exceeded  its  range  of  applicability. 


Test  #3  Simulation,  v=2147  m/s  @  15°. 

Figure  87  through  Figure  89  show  the  simulation  at  interesting  times  of  2.5,  5.5, 
and  10.0  microseconds. 
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Figure  87  -  Test  #3  Sim,  t=2.5  jus  Figure  88  -  Test  #3  Sim,  t=5.5  jus 
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Figure  89  -  Test  #3  Sim,  t=10.0  jus 
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As  with  Test  #2,  from  2.5  to  5.5  ps,  the  projectile  and  target  start  do  deform  and  form  a 
gouge.  At  10  ps  the  projectile  has  formed  two  cracks,  one  on  the  top  surface,  and  one  on 
the  bottom  surface.  The  penetration  depth  of  the  target  is  roughly  1  cm.  To  determine  a 
better  estimate,  Figure  90  shows  a  blown  up  image  of  the  simulation  at  10  ps.  Note  that 
the  x  and  y  axis  scaling  are  different. 
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Figure  90-  Test  #3  Sim,  t=10.0  /us 

It  is  seen  that  the  simulated  penetration  depth  is  on  the  order  of  1.15  cm.  The  measured 
depth  is  1.0  ±  0.1  mm.  The  difference  in  depth  is  roughly  13  percent.  Pressure  and 
plastic  strain-rate  profiles  are  similar  to  Test  #2  and  are  shown  in  Appendix  C. 

Bringing  this  altogether,  Table  12  shows  the  comparison  of  the  measured  and 
simulated  depths  of  the  experiments. 
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Table  12  -  Comparison  of  Experimental  and  CTH 


Test 

Velocity 

(m/s) 

Angle 

(deg) 

Uy 

(m/s) 

Experiment 

CTH  Simulation 

Absolute 

Difference 

(mm) 

Percent 

Difference 

Simulated  Depth 

(mm) 

1 

2225 

■a 

386 

0.5 

0.7 

0.2 

29% 

n 

2150 

556 

1.0 

1.25 

0.25 

20% 

2147 

■9 

556 

1.0 

1.15 

0.15 

13% 

n 

2163 

10 

376 

0.5 

0.7 

0.2 

29% 

Results 


The  results  show  that  while  CTH  cannot  fully  capture  the  entire  event  in  time,  the  current 


model  does  yield  accurate  results.  With  strain  rates  up  to  104  /s,  the  use  of  the  current 


strength  model  has  exceeded  its  range  of  applicability. 


Equivalent  HHSTT  Sled  System 

Using  the  methodology  described  earlier,  the  equivalent  HHSTT  can  be 
developed  using  the  Rod  Impact  results  from  Table  1.  Test  #1  is  used  to  give  the 
equivalent  HHSTT  sled  system.  It  should  be  noted  that  only  the  velocity  and  angle  of  the 
equivalent  system  are  affected  by  the  different  tests.  Table  13  compares  the  two  models 
side-by-side. 


Table  13  -  Model  Comparison 


Parameter 

Rod  Impact 
Model 
(Test  #1) 

Equivalent 
HHSTT  Sled 
System 

m  (g) 

4.78 

2.00E+05 

1  (cm) 

2.5 

20.32 

h  (cm) 

0.55 

4.47 

w  (cm) 

- 

10.80 

v  (m/s) 

2225 

812.35 

ux(m/s) 

2192.2 

800.00 

uy  (m/s) 

-386.37 

-141.10 

It  is  worth  mentioning  here  that  this  equivalent  system  doesn’t  exhibit  typical 
HHSTT  scenarios.  The  nominal  HHSTT  slipper  height  is  2.5  cm,  whereas  the  equivalent 
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system  height  is  4.47  cm.  This  is  due  to  the  fact  that  the  aspect  ratio  is  constrained  by  the 
Buckingham  Pi  Theorem.  Also,  the  nominal  impact  angle  is  roughly  0.03°,  whereas  the 
equivalent  system  impacts  at  10°.  Again,  this  is  due  to  fact  that  the  impact  angle  must 
remain  the  same  between  models. 

Using  the  time  scale  of  Equation  65,  every  1  time  unit  of  the  Rod  Impact  model 
equals  22.3  time  units  of  the  equivalent  HHSTT  sled  model.  With  the  simulation  of  Test 
#1  running  to  10  ps,  the  equivalent  HHSTT  sled  model  should  be  simulated  out  to  223 
ps.  Unfortunately,  modeling  this  larger  system  out  this  far  requires  much  more 
processing  power  than  currently  available  here  at  AFIT.  Because  of  this,  the  equivalent 
model  is  run  out  to  only  10  ps.  Figure  91  through  Figure  93  show  the  displacement, 
pressure,  and  strain-rate  profiles  at  10  ps. 


Material  impact 


0  5  10  15  20  25 

2DR  X  (cm) 

Equilalwit  HHSTT:  *=812.35,  anq1e=101  ux=B00,  uy=  141.1  m/s 


0  5  10  15 

2DR  X  (cm) 

ttquNcIgnf  HHSTT:  v=B12,35(  «igEe=K31  ux=800,  uy= 


Figure  91  -  Equivalent  HHSTT  Model,  t=10  jus 


Figure  92  -  Pressure,  Equivalent  Model,  t=10  jus 
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Psnatic  Strain  Rate 


He  Strain  Rate 


ujowd1 

fl.ElJit)1 

7K&9*HJ-’ 

7.MM03 

G.KMjcW* 

e.ii5*K>J 

5.727<YH 

.MWV*)3 

-UtiKW3 

CPW 

AJMftrtO* 

iHJiKW5 

iSIStB3 

5-iSnpiO1 

1JB12<I0'5 

1.454.Y1&1 

6.7GSmHP 

iSShtt* 

-LODD.B1 


0  5  10  15  20  25 

20  R  X  (cm) 

FquSIfjlant  HHSTT:  v=8l2.35a  angt^tt),  ux=S0O,  uy=-141,1  m/s 


Figure  93  -  Strain-Rate ,  Equivalent  Model,  t=10  jus 


As  can  be  seen,  pressures  of  3  GPa  and  strain-rates  up  to  104  /s  are  present  just  as  in  the 


Cinnamon  simulations.  Figure  94  shows  a  close-up  of  the  impact  zone,  note  the  y-axis 


scaling. 
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■  wi 
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10  12  14  16  18  ZO  22 

ZDR  X  (cm) 

Eqpilnlent  HHSTT:  v=812.35,  angle  =1Ga  ux=800„  uy=  141.1  m/s 

Figure  94  -  Equivalent  HHSTT  Model,  t=10  jus 


As  can  be  seen,  the  simulation  the  equivalent  model  predicts  a  penetration  depth  of 


roughly  1.0  mm  at  10  ps.  Although  the  CTH  model  cannot  predict  penetration  depth 


beyond  10  ps,  Figure  94  does  show  the  initial  signs  of  gouge  formation  through  the  high 
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pressure  core  and  material  interaction  at  the  lower  front  area  of  the  slipper,  suggesting 
that  the  equivalent  HHSTT  sled  system  will  eventually  develop  a  classic  gouge. 

This  suggests  that  although  the  test  facilities  here  at  AFIT  are  not  capable  of 
velocities  in  the  realm  of  the  HHSTT,  it  is  possible  to  utilize  these  facilities  to  further  the 
understanding  of  this  gouging  phenomenon.  With  the  Buckingham  Pi  and  CTH  models 
producing  good  results,  it  is  possible  to  relate  the  experiments  made  here  to  real  sled 
systems  at  the  HHSTT. 


V.  Conclusions  and  Recommendations 


Conclusions 

The  primary  goal  of  this  research  was  to  compare  a  computer  simulation  with 
experiment  of  an  oblique  ballistic  impact.  Through  the  Buckingham  Pi  Theorem,  this 
model  was  scaled  up  to  an  equivalent  HHSTT  sled  system.  While  this  equivalent  system 
does  not  physically  exist,  it  did  give  ballpark  estimates  of  what  similar  systems  should 
experience  under  the  same  conditions  as  the  experiments.  The  hydrocode  computer 
program  CTH  was  used  to  validate  that  it  can  simulate  a  cylindrical  projectile  impacting  a 
target  at  an  oblique  angle  using  a  plane-strain  analysis;  verify  the  experimental  ballistic 
impacts;  and  to  simulate  the  equivalent  HHSTT  sled  model  to  see  if  the  gouging 
phenomena  were  present.  To  reach  the  primary  goal  of  this  work,  a  series  of  model 
comparisons  were  made.  The  conclusions  of  this  work  are  summarized  in  the  following 
bullets. 

Through  comparisons  of  CTH  simulations  of  axisymmetric  and  plane-strain 
vertical  impacts,  it  was  shown  that  CTH  can  accurately  simulate  an  axisymmetric 
impact  event  using  a  plane-strain  analysis  up  to  10  ps.  Pressure  profiles 
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comparisons  showed  that  they  differ,  but  only  by  15  to  30  percent  at  most  in 
magnitude. 

It  was  shown  that  the  predicted  penetration  depths  are  greater  than  the  Cinnamon 
experiments  by  13  to  29  percent,  which  is  quite  excellent.  Pressures  and  strain 
rates  of  3  GPa  and  104  /s  were  observed.  This  shows  that  while  the  Johnson-Cook 

'J 

strength  model  was  based  on  strain  rates  up  to  10  /s,  the  current  model  is 
accurate. 

Through  the  use  of  dimensional  analysis,  an  equivalent  HHSTT  sled  system  was 
developed  based  on  Test  #1  experiment.  Unfortunately,  modeling  this  larger 
system  out  to  the  required  time  required  much  more  processing  power  than 
currently  available.  However,  the  simulation  was  run  out  to  10  ps,  and  showed 
signs  of  gouge  initiation  as  in  the  early  stages  of  the  Cinnamon  simulations. 

This  suggests  that  although  the  test  facilities  here  at  AFIT  are  not  capable  of 
velocities  in  the  realm  of  the  HHSTT,  it  is  possible  to  utilize  these  facilities  to 
further  the  understanding  of  the  gouging  phenomenon.  Since  the  Buckingham  Pi 
and  CTH  models  produced  good  results,  it  is  possible  to  relate  the  experiments 
made  here  to  real  sled  systems  at  the  HHSTT. 

Recommendations  for  Future  Research 

This  investigation  assumed,  along  with  Szmerekovsky,  Rickerd,  and  others,  that 
the  mass  of  the  HHSTT  sled  system  is  evenly  distributed  over  the  four  slippers,  and  that 
the  slipper  is  oriented  flat  and  level  against  the  rail.  While  this  may  be  a  fairly  good 
assumption  when  the  sled  is  at  rest,  it  is  most  likely  not  true  at  high  velocities. 
Aerodynamic  forces,  rail  characteristics,  and  sled  dynamics  can  have  profound  effects  on 
weight  distribution  and  pitch  angle.  The  validity  of  this  assumption  should  be 
investigated. 

The  HHSTT  has  an  epoxy  coating  on  the  rails,  which  serves  as  a  type  of  lubricant 
to  mitigate  gouging.  Neither  metallurgical  analysis  on  the  experimental  target  rails  nor 
effects  of  coatings  was  within  the  scope  of  this  work.  Addition  of  coatings  to  the  CTH 
models  would  allow  improved  comparisons  to  the  actual  experiments.  This,  along  with  a 
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metallurgical  analysis  of  the  target  rails,  would  enhance  the  understanding  of  the  HHSTT 
gouging  phenomena. 

The  Johnson-Cook  coefficients  used  in  this  work  was  based  on  the  prior  work  of 
Kennan  and  Yun.  They  utilized  Split  Hopkinson  Bar  tests  in  their  analysis,  which  were 
only  able  to  give  data  on  strain  rates  up  to  103  /s.  However,  strain  rates  on  the  order  of 
104  /s  were  seen  in  the  CTH  simulations.  To  more  accurately  define  the  strength  models, 
other  testing  methods  that  are  capable  of  yielding  data  on  higher  orders  of  magnitude  of 
strain  rates  should  be  utilized,  such  as  Flyer  Plate  tests. 

The  impact  angles  in  the  Cinnamon  experiments  were  chosen  two  to  three  times 
the  minimal  angle  required  to  theoretically  gouge  the  rail.  This  was  done  to  ensure 
gouging  of  the  rail.  However,  it  was  shown  that  the  higher  the  vertical  velocity,  the 
shorter  the  time  frame  CTH  can  simulate  a  cylindrical  projectile  impacting  a  target  at  an 
oblique  angle  using  a  plane-strain  analysis.  This  resulted  in  CTH  not  being  able  to  fully 
simulate  the  event  in  time.  This  also  led  to  an  equivalent  HHSTT  model  with  the  same 
impact  angle,  which  is  not  typically  seen  in  the  field.  More  impact  experiments  should 
be  utilized  with  a  compromise  between  matching  impact  angle  and  vertical  kinetic 
energy,  since  it  is  difficult  to  match  both  at  the  same  time  here  at  AFRL. 
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Appendix  A.  CTH  Input  Deck  of  Axisymmetric  -  Plane  Strain  Comparison 


*eor*  genin 

Axisymmetric  Impact:  V300  on  1080,  v=375  m/s 

control 

mmp 

ep 

vpsave 

endcontrol 


mesh 

block  1  geom=2dc  type=e 


X 

o 

II 

o 

0 

xl 

n=60 

w=0 .6  dxf=0 . 01 

x2 

n=8 

w=0 . 4  dxf =0 . 05 

x3 

n=2  0 

w=2 .  dxf=0 . 1 

endx 

i 

ll 

o 

>1 

.0 

yi 

n=2  0 

w=2  dyf=0 . 1 

y2 

n=2  0 

w=l  dyf =0 . 05 

y3 

n=350 

w=3 . 5  dyf =0 . 0 1 

endy 

endb 

endmesh 


2dc  is  two  dimensional  cylindrical 
e  is  for  an  Eulerian  solution 


insertion  of  material 


block  1 

package  topblock 
material  1 
numsub  50 
yvel  -375.e2 

insert  box 


* 

* 

* 

* 

* 

* 

* 


* 

pi  0.0  0.275  * 

p2  0.275  2.5 
endinsert 


Change  only  the  first  number,  leave  'e2' 
this  converts  m/s  to  cm/ s . 

This  is  where  you  input  the  cylinder  - 
only  model  1/2  of  cylinder,  and  then 
"mirror"  in  plotting. 

Format  is  pi  is  the  bottom  center  point  of 
cylinder . 

p2  is  top  right  hand  corner  of  cylinder. 
Change  pi  and  p2  to  define  the  size  of  rod. 


delete  circle 
center  0  0.275 
radius  0.275 
enddelete 
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endpackage 


package  tip  *  Insert  circular  rod  tip 

material  1 
numsub  50 

yvel  -375.e2 

insert  circle 
center  0  0.275 
radius  0.275 
endinsert 

endpackage 

package  target  *  Insert  Target 

material  2 
numsub  50 
insert  box 
pi  0.0,  0.0 

p2  3,  -4.0 
endinsert 
endpackage 
endblock 
endinsertion 

edit 
blockl 
expanded 
endblock 
endedit 

tracer  *  Tracer  Points 

add  0  0.01 

add  0  0.275 

add  0.275  0.01 
add  0.275  0.275 
endt 

eos 

*  Information  for  metals 
MAT1  SES  STEEL_V300 

CTH 

MAT 2  SES  IRON 
enough  for  1080 
endeos 

epdata 
mix  3 

matep  1  JO  USER  *  J-C  coefficients  for  VascoMax  300. 

ajo=  2.17el0  *  Dynes/cmA2 

bjo=  0.124el0  *  Dynes/cmA2 

c j  o=  0.004  6 
mj  o=  0 . 95 
njo=  0.3737 
tjo=  . 040161e-l 


*  EOS  for  Vascomax  300  already  in 

*  EOS  for  Iron  should  be  close 
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poisson=  0.283 

matep  2  JO  USER 
ajo=  5.25e9 
bjo=  3.59e7 
c j  o=  0.29 
mjo=  6.525e-l 
n j  o=  0.6677 
tjo=  . 1591885e-l 
poisson=  0.27 

vpsave 
lstrain 
endep 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

*eor*  cthin 

Axisymmetric  Impact:  V300  on  1080,  v=375  m/s 
control 
mmp 

tstop  =  15.0e-6  *  you  may  need  to  increase  the  stop  time  to  reach 

*  end  of  event. 

endc 


*  J-C  coefficients  for  1080  steel. 

*  Dynes/cmA2 

*  Dynes/cmA2 


Convct 

convection=l 

inter f ace=high_resolution 
endc 


-2 . 4el0 
-1 . 38el0 
-12 . 0e9 
-12 . 0e9 


f racts 
pf racl 
pf rac2 
pfmix 
pfvoid 
endf 

edit 
shortt 
tim  0.0, 
ends 
longt 

tim  0.0, 
endl 
plott 
tim  0 . 0 

endp 
histt 
tim  0.0, 

htracer  all 
endh 
ende 


dt  =  1.0 

dt  =  1.0 

dt  =  0.25e-6 

dt  =  0.25e-6 


*  this  sets  the  time  step  it  plots  at 

*  1/2  microsecond 


*  this  sets  the  time  step  it  plots  at 

*  1/2  microsecond 
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boundary 

bhy  *  Rigid  boundaries  all  around 

bl  1 

bxb  =  0  ,  bxt  =  1 
byb  =  1  ,  byt  =  1 
endb 
endh 
endb 

cellthermo 

mmp3  *  This  was  recommended  by  Eglin  and  appears 

*  to  give  good  results  as  well. 

ntbad  1000000 
endc 
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Appendix  B.  CTH  Input  Deck  for  Rod  Impact  Test  #1 


*eor*  genin 

Oblique  Impact:  V300/1080,  v=2225,  angle=10,  ux=2192.2,  uy=-386.37  m/s 

*  All  output  from  CTHgen  will  have  this  title. 

* 

*  This  is  an  impact  scenario  of  a  Vascomax  300  projectile 

*  hitting  a  1080  steel  target. 

* 

control 

mmp 

ep 

vpsave 

endcontrol 

* 

* 

* 

mesh 


block  1 

geom 

=2dr 

type=e 

* 

This  section  defines  the  mesh,  not 

* 

the  cylinder  or  block. 

* 

2dr  =  2D  rectangle  (plain  strain) 

* 

Use  2dc  for  axisymmetric . 

* 

e  =  Eulerian  solution. 

X 

o 

II 

o 

0 

* 

xO  =  left  starting  value  of  x. 

xl 

3 

ii 

i—1 

o 

o 

i — i 

ll 

£ 

dxf=0 .100 

* 

xl,x2  =  x  subzone  for  defining 

* 

different  meshes  along  x. 

x2 

n=50 

LO 

o 

II 

£ 

dxf=0 . 010 

* 

ALL  DIMENSIONS  MUST  BE  IN  cm! 

x3 

n=22  0 

\ — 1 

\ — 1 

II 

£ 

dxf=0 . 005 

* 

n  =  #  of  cells;  w  =  section  width; 

x4 

n=4  0 

w=0 . 4 

dxf=0 .010 

* 

dxf  =  dx  first;  dxl  =  dx  last 

* 

(optional) 

x5 

n=2  0 

o 

C\] 

II 

£ 

dxf=0 . 100 

endx 

i 

ll 

o 

>i 

.0 

* 

y  has  same  inputs  as  x 

yi 

n=2  0 

o 

CN] 

II 

£ 

o 

\ — i 

o 

II 

4-1 

>i 

T3 

* 

yO  =  bottom  starting  position. 

y2 

n=140 

w=l .  4 

dyf =0 . 0 1 

y3 

n=2  4  0 

w=l .  2 

dyf=0 .005 

y4 

n=4  0 

w=0 . 4 

o 

\ — i 

o 

II 

4-1 

>i 

T3 

endy 

endb 

endmesh 

* 

* 

* 

insertion  of  material 
block  1 

package  rod  *  Insert  rectangular  rod  and  rotate 

material  1 
numsub  100 


80 


xvel  2192. 2e2  *  Change  only  the  first  number,  leave  'e2' 

*  this  converts  m/s  to  cm/ s . 

yvel  -386.37e2 

insert  box 

pi  0.0  0.0 
p2  2.225  0.55 


endinsert 

delete  circle  *  Delete  a  circle  for  tip  insertion  (next  step) . 
center  2.225  0.275 
radius  0.275 
enddelete 
endpackage 

package  tip 
material  1 
numsub  100 
xvel  2192 . 2e2 
yvel  -386.37e2 

insert  circle 

center  2.225  0.275 
radius  0.275 
endinsert 
endpackage 

package  target  *  Insert  Rail  Material 

material  2 
numsub  100 
insert  box 
pi  0.0  -4.0 
p2  5.0  0.0 
endinsert 
endpackage 

endblock 
endinsertion 
* 

* 

* 

edit 
blockl 
expanded 
endblock 

endedit 
* 

* 

* 

eos 


*  Eg  Of  State 


*  Insert  circular  rod  tip 


*  Separate  package.  Must  include  velocity. 


*  This  is  where  you  input  the  projectile. 

*  Format:  pi  =  bottom  center  point  of 

*  cylinder. 

*  p2  =  top  right  hand  corner  of  cylinder. 

*  Change  pi  and  p2  to  define  the  size  of  rod. 
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MAT1  SES  STEEL_V300 

MAT 2  SES  IRON 
endeos 


*  Sesame  Eq  of  State  tables,  limited 

*  materialselection . 

*  Iron  =  closest  material  available  to  1080. 


epdata 
mix  3 


*  Elastic/Plastic  data 

*  mix  3  =  normalized  vol  avg  yeild  strength 


matep  1  JO  USER 
ajo=  2.17el0 
b j  o=  0 . 124el0 
c j  o=  0.004  6 
mj  o=  0 . 95 
n j  o=  0.37  37 
t j  o—  . 040161e-l 
poisson=  0.283 


*  J-C  coefficients  for  VascoMax  300. 

*  Dynes/cmA2 

*  Dynes/cmA2 


matep  2  JO  USER  *  J-C  coefficients  for  1080  steel. 

ajo=  5.25e9  *  Dynes/cmA2 

bjo=  3.59e7  *  Dynes/cmA2 

c j  o=  0.29 
mjo=  6.525e-l 
n j  o=  0.6677 
t j  o—  . 1591885e-l 
poisson=  0.27 


vpsave 

lstrain 

endep 

* 

* 

* 

tracer 

*  Tracer  input,  starting  at  xl,yl  to  ending  x2,y2. 

*  n=number  of  tracers  to  distribute  including 

endpoints . 

*  Projectile  boundary 

add  1.0  .01  to  3.0  .01  n=50 

endt 

* 

* 

* 

*eor*  cthin 


Oblique  Impact:  V300/1080,  v=2225,  angle=10,  ux=2192.2,  uy=-386.37  m/s 


control 

mmp 

tstop  =  15.0e-6 
cpshif t=600 . 
endc 
* 

* 


*  Defines  sys  parameters  for  execution. 

*  mmp  =  multiple  pressures  and  temps  in  mixed 

*  material  cells. 

*  tstop  =  Problem  stop  time  in  seconds . 

*  cpshif t  allows  extra  time  for  CTH  to  right  data. 
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* 

Convct 

convection=l 


inter f ace=high_resolution 

endc 

* 

* 

* 


*  Convection  control  input. 

*  Convect  internal  energy  based  on 

*  energy  density  and  mass  density. 

*  Discards  kinetic  energy. 

*  interface  tracker,  high_res 

*  recommended  for  2D. 


f racts 

pfracl=-2 . 4el0 
pfrac2=-1.38el0 
pfmix=-12 . 0e9 

pfvoid=-12 . 0e9 

endf 

* 

* 


*  Fracture  Strengh 

*  pfrac#  =  fracture  pressure  or  stress  of  material  # 

*  pfracl  and  2  verified  by  J  Cinnamon. 

*  pfmix  =  fracture  stress  or  pressure  in  cell  with 

*  mixed  mat  no  voids . 

*  pfvoid  =  fracture  press  or  stress  in  cell  with 

*  void. 


* 

edit 

shortt 

tim  0.0,  dt=l . 0  *  Restart  data  will  be  written  every  ' dt ' 

*  seconds . 

ends 


longt 

tim  0.0,  dt=l . 0  *  Restart  data  will  be  written  every  1 dt * 

endl 


plott 

tim  0.0  dt=0.25e-6 
endp 


*  cthplot  data  is  written  to  restart  file 

*  every  ' dt '  seconds  starting  at  'tim'. 

*  Beware  -  Restart  file  size  limited  to  2GB. 


histt 

tim  0.0  dt=0.25e-6 
htracer  all 

endh 


*  History  data  will  be  written  to  hcta  every 

*  1 dt *  seconds.  Data  will  be  written  for  all 

*  tracer  points . 


ende 

* 

* 

* 

boundary 

bhy  *  rigid  boundaries  all  around 

bl  1 

bxb  =  1,  bxt  =  1 
byb  =  1,  byt  =  1 
endb 
endh 
endb 
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* 

* 

* 

vadd  *  Apply  vel  to  keep  gouge  in  mesh 

block=l 
tadd=0 . 0 
xvel=-2 1 92 . 2e2 

endvadd 
* 

* 

* 

cellthermo 
mmp3 

ntbad  1000000 
endc 
* 

* 

* 


*  This  was  recommended  by  Eglin  and  appears 

*  to  give  good  results  as  well. 
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Appendix  C.  CTH  Plots  of  Cinnamon  Simulations 


Test  #4: 


Pressure 


Plastic  Strair>-Rale 


Ptoiit  itrah  km* 
(/*) 


Oblique  Impact:  V3G0/10BQ,  v=2163,  angl&=10,  ux=2  80.14,  uy  =-375.6  m/s 
BMCEFN  2/13/06  16:5686  CTH  3142  Tlme=l. 00003x10“ 5 


0 . 0 
20R 

Oblique  Impact:  V3G0/1DB0,  v=2163s 
BMCEFN  2/13/OS  13:56:36  CTH 


=2130.14,  uy=-375,E 
Time  “  1.0000  3 >  10“ 


Test  #3: 


Pins  tic  Stroin— Rate 


0  .  O 

2DR  X  (cm) 

l  Impncl;  V3OP/10BQ,  v=?147,  anglB=15,  0*^073.34, 


2/13/06  16:57:55  CTH 


3161  Tims -1.00031x10' 


535-68  m/g  Qb |:qLJfl  Impcigl:  V3OD/10BO, 


2/li/OE  IE :S?:£5  C1H 


-555.68  m/5 


31E4  Tls™  =  l.O0Ci!x1O^ 
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